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Summary 


The possibility of predicting the ultimate static bearing 
capacity of a pile from its dynamic behavior during driving is 
investigated by the wave equation method of pile driving analysis. 
The method developed by E. A. L. Smith over a period of years is 
followed, but its validity is first shown by proving that the math- 
ematical model used is equivalent to the wave equation. Computer 
programs developed for purposes of this investigation are presented 
in their entirety together with detailed instructions on their use 
so that they may be utilized by future investigators who may not be 
familiar with computer work. 

The soil engineering aspects of the problem are explored. The 
information available on dynamic soil properties from published test 
data is reviewed to obtain information applicable to pile driving. 
Using values thus deduced and the computer programs, computations of 
ultimate resistance are made from 2) published pile driving records 
and correlated with their load test results. 

Although not conclusive because of the relatively small num- 
ber of correlations attempted, the results are very encouraging. It 
appears that the wave equation method of pile driving analysis may 
become an accurate tool in predicting a pile's static bearing capacity 
from its driving record. This appears to be the case regardless of 
the type or size of pile and driving equipment, and indications are 
that predictions may be possible for piles driven in cohesive as well 


as granular soils. General limitations of the method are summarized. 
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Values deduced for point and side damping and ground quake 
in various soils are presented. Corresponding values of friction 
acting on the pile sides as a percentage of the ultimate resistance 
are given for the piles investigated. More research into this problem 
is indicated, and recommendations are made for such work, 

It is concluded that even with inexact knowledge of dynamic 
Soll properties the wave equation method of pile driving analysis 


appears to give good results. 
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Chapter I 


INTRODUCTION 


The possibility of relating the behavior of a pile during 
driving to its static load carrying capacity by the use of the wave 
equation method of pile driving analysis is encouraging. Many 
engineers today realize that the numerous dynamic pile driving 
formulas which have been used, and are still being used, have seri- 
ous limitations and cannot be depended upon to give reliable results, 
but they are still used for lack of an adequate substitute. Some 
feel that the dynamic formula approach is not valid, and have turned 
instead to trying to determine the bearing capacity of a pile from 
purely static soil mechanics and structural considerations for the 
cases where the pile-soil system is required to carry only static 
loads. Such an approach is certainly valid, and it requires a thorough 
investigation of the soil conditions, together with a careful evalua- 
tion of their engineering properties, and the use of a great deal of 
personal judgment in its application. 

It is believed that the wave equation method of analysis of 
the dynamic pile driving problem is an equally valid approach and may 
offer advantages over the static approach if it can be successfully 
applied. This is not to say that one approach should be used to the 
exclusion of another, but that they are working toward a common goal 
and should supplement each other, The wave equation method is not 
just another dynamic pile driving formula, but, rather, itis a 


method of analysis which is well founded mathematically. It requires 
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a knowledge of the soil and its properties, both static and dynamic, 
for its successful application. Current knowledge of soil dynamics 
is incomplete, but an increasing amount of work is being done in this 
preld, 

Approaching the problem from the dynamic side, the over-all 
consideration is to relate the dynamic behavior of the driving equip- 
ment-=pile-soil system to the static behavior of the structure-pile- 
soil system with due consideration to changes in properties of the 
soil during and after driving of the piles. A part of this broader 
problem is to relate the dynamic resistance to driving of a single 
pile under the last hammer blow to its static bearing capacity, with 
time effects between driving and static loading minimized. It is for 
this more limited problem that the wave equation method of analysis 
aceuserul. 

The wave equation alone will not provide a solution to the more 
general consideration which includes the group effect of piles, de- 
velopment of negative friction, consolidation of a clay layer beneath 
the pile tip, long-term changes in water table, deterioration of the 
pile due to insect attack or deleterious chemicals, set up or relaxa- 
tion of the soil with time, or development of hydrostatic uplift. 
That it does not provide a solution to such problems may be obvious, 
but the statement is made here in an effort to clear up the type of 
misunderstandings and objections already voiced in opposition to the 
application of the wave equation to pile driving analysis. In short, 


it is not a magical formula, but simply a tool which is well founded 
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mathematically and is available to assist in visualizing and evaluating 
the dynamics of the problem. 

The wave equation can be utilized to advantage for investigation 
of other facets of pile driving such as stresses in the pile and 
selection of appropriate pile driving equipment for a particular set 
of field conditions; but the main points of this investigation will 
be to consolidate the wave equation theory as applicable to pile driv- 
ing, make readily available computer programs to facilitate the work 
of future investigators, explore the unsolved problems of interaction 
between soil and pile, attempt correlation between wave equation solu- 
tions and a few pile driving records and load tests, and to suggest 


possibilities for further research work in this general area, 





Chapter II 


BACKGROUND 


Although use has been made of piles for at least 2,000 years, 
attempts to determine the bearing capacity of a pile from its driving 
record have been generally unsuccessful in spite of the amount of 
effort spent in this direction. An extensive treatment of the develop- 
ment and interrelatedness of dynamic formulas by Chellis (ref. 5) is 
included, and he shows that they may be grouped as follows: 

The first dynamic formula was proposed by Major Saunders in 
1851 and is perhaps the simplest and most direct. He equated the 
weignt of the ram multiplied by the stroke,with the driving resistance 
multiplied by the set, and applied a factor of safety of 8. 

Another group of proposed formulas contained a fixed coefficient 
which was supposed to compensate in some degree for factors affecting 
the results which were not included as terms in the formulas. The 
Engineering News, Wellington, Vulcan, and Bureau of Yards and Docks 
(not now used by Navy) formulas are of this type. 

Another group of formulas attempts to account for the variables 
by using expressions for efficiency of applied energy by including 
relative weights of the pile and the hammer. The Dutch, Ritter, and 
Benabencq formulas are of this type. 

Other similar formulas try to include the effect of variables 
by using both fixed coefficients and expressions for the relative 
weights of the pile and ram. The Eytelwein and Navy-McKay (not now 


used by the Navy) are of this type. 
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The next group contains either all or part of a series of 
terms designed to represent impact losses in the driving cap, soil, 
and pile during driving. The Redtenbacher and Hiley formulas are 
among those which fall into this category. 

Only a few of the formulas proposed have been mentioned. 
Chellis (ref. 5) lists 38, while it is understood that the Editors of 
Engineering News Record have on file 450 such formulas. Some of these 
formulas are extremely simple, while others are quite complicated. 
Some try to approximate the dynamic situation, while others are con- 
trived from a purely statistical approach without reference to the 
numerous variables actually involved. For example, Marvin Gates has 
proposed (ref. 1) that for a pile driven with a Vulcan Number 1 Hammer 
the ultimate bearing capacity may be determined by multiplying 8 by 
the log of (10/set in inches). He claims, ". . .this relationship 
gives more consistently accurate results than the most complex dynamic 
formula yet advanced," in spite of the fact that it contains but two 
simple parameters. That this can be true for a formula which seems 
to ignore the basic aspects of the problem is not surprising when one 
considers the shortcomings and omissions of the other dynamic formulas 
proposed. 

Some of the dynamic formulas were based on Newtonian theory 
of impact which has been shown by A. E. Commings (ref. 10) as in- 
applicable to the pile driving problem. Other formulas subtracted the 
same energy losses twice,as also shown by Cummings (ref. 10). The 
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Society of Civil Engineers recognized the many limitations of dynamic 
pile driving formulas (ref. 8). After this point of agreement, en- 
gineers seem to be divided into two general schools of thought. One 
group is optimistic and feels that rational use can be made of such 
formulas. The other group stresses the limitations of the formulas, 
and basically takes the position that, since none of them are reliable, 
you may as well save yourself work and choose a simple formula. The 
Jatter school of thought may seem rash at first reading; but it does 
have some merit, for it is simply recognizing that all of the dynamic 
formulas are over-simplifications of a very complex problem. 

It is certain that each of the 50 formas available to 
choose from will give reliable results under particular combinations 
of driving equipment and soil conditions for a particular pile, but 
there would seem to be an almost infinite number of combinations of 
the variables involved making the choice of a formula a difficult one. 
That so many formulas have been developed is indicative of the interest 
in the problem and the ingenuity of engineers in attempting to estab- 
lish a relationship between static and dynamic resistance of piles. 
Empirical formulas certainly have their place when their range of 
applicability is clearly established and when their limitations are 
kept clearly in mind; however, this has not always been the case in 
pile driving. 

In spite of the interest exhibited, there has not been a great 
deal of work done theoretically or experimentally in defining or 


attacking the fundamental problems involved. The complexity of the 
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problem, the difficulty in isolating the variables, and the limited 
techniques and equipment available for detailed analysis have generally 
limited advances in this field. Full utilization, however, of the 
facilities available has not been realized. 

The basic problem is to relate the dynamic resistance to 
driving of a single pile to its static bearing capacity with time 
effects between driving and static loading minimized. Available to 
assist in the analysis of this problem is the wave equation which was 
developed over a hundred years ago by De Saint Venant and Boussinesq 
for end impact on rods. It is well founded mathematically. D. V. 
Issaccs, in 1931, was the first to point out that wave action occurred 
during the driving of piles. In 1938, E. N. Fox published a solution 
to the wave equation applied to pile driving, but, as no electronic 
computers were available at that time, he made a number of simplifying 
assumptions. Again in 190 and 191, A. E. Cummings reported on the 
work of the foregoing investigators. Ten years elapsed before the wave 
equation was again proposed as being applicable to pile driving analysis, 
this time by E. A. L. Smith (ref. 27). Then in 1961 Smith published a 
more comprehensive article on application of the wave equation to pile 
driving and described a numerical approach together with a description 
of his efforts in developing a computer program for the solution (ref. 
28). His efforts have been untiring in developing this approach and 
have opened the door for much needed work in this field, 

The wave equation, as will be shown, mimics the behavior of 


a pile and provides a good indication of what will happen to the pile 
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under conditions of impact for specific boundary conditions. For this 
reason, a distinction is made between the wave equation and dynamic 
formulas which appear to be mathematical but are in reality empirical. 
The present state of knowledge up through Smith's work can provide 

a prediction of effects caused by a ram hitting a pile. Still to be 
investigated 1s the nature of the effects caused by interaction be- 
tween the pile and the soil. This, too, is a dynamic problem, and 

it lies in the new field of soil dynamics. This matter will be ex- 
plored, but, first, the wave equation theory and its numerical solution 


will be developed. 
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Chapter III 


DEVELOPMENT OF BASIC EQUATIONS AND THEORY 


A, Theory 


A. E, L. Smith, over a period of severalyears, derived equa- 
tions suitable for numerical calculations and computer adaptation 
for pile driving analysis by the wave equation. His derivations are 
based on a mathematical model using weights and springs. This work 
has generally been accepted as correct and as equivalent to the wave 
equation,,. but the full derivation establishing this equivalence has 
not been published. In order to establish this equivalence and the 
validity of the numerical method, a step-by-step derivation from first 
principles is presented. For completeness, Smith's derivations from 
the mathematical model are included, 

The general plan for establishing this link is as follows: 

a. <A partial differential equation is developed and is shown 
to be a form of the wave equation. 

b. This equation is converted into a difference equation for 
ease of numerical solution. 

c. To simplify further numerical solution, a mathematical 
model consisting of weights and springs is assumed as being equivalent 
to the elastic rod used as a basis for steps a and b above. Five 
equations are developed from this mathematical model which is based 
on the assumption that all springs are perfectly elastic,and the pile 


is represented typically as shown by figure (3). 













— eine (aos ae as 
— io ee = am = 
ee ey Srakcangilttadl hee 
Se rte 0 el ol fee ot 
nA ent Cee Fre ye ey AD Re hb lle - 
re Eee elt eae il et 
yah IS eg, Cintra 
ee ee | 
Nee fe tes ye IE Ant ow 
hee” meee ae tee b derermesmeceeme + prorrrL A emm am 6h 
100 ee ha ae emma, 5 Camm Samm os a 
—— (8! « pateeéel ace te Ore yale ape 
-t —- weer’ ~s] eons — een 

4) ay Gli 6 4) Oxi eel tue eemarearel 
1.) ©) ey Pewee = ee ~ Ug Ft ee 
(| aes 6 Cee tf Genki eerie 











10 


d. It is then shown that the equations of step c.can be com- 
bined to form the same equation developed in step b. above, and it is 
concluded that the mathematical model and the five equations developed 
from it are equivalent to the elastic rod, and the wave equation is 
applicable to it. 

Beginning with development of the partial differential equation, 
as shown by Timoshenko and Goodier (ref, 32), consider the longitudinal 
waves in a prismatic bar of constant cross section such as shown in 
figure (1). The axis of the bar is taken as the x-axis. It is assumed 
that: 

a. Cross sections of the bar remain plane during deformation. 

b, The unit elongation at any cross section mn, due to a 
longitudinal displacement D, is equal to . and the corresponding 
tensile force in the bar is AE 3e , where A is the cross sectional 
area and E is Young's modulus. 

c, A simple tension in the x direction and the unit elongation 


= is accompanied by lateral contraction of the amount ~v = where -y 


ax? 
is poisson's ratio. Inertia forces corresponding to lateral particle 
motion are neglected since the length of the waves is large in comparison 
with the cross sectional dimensions of the bar. 


The foregoing leads to the equation for force at any section 


in the bar of uniform cross section: 


\ 
P=FaAcoD, 
ag X 
The net force on a section of bar dX, figure (1), is: 
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assuming that E may not be constant for purposes to be shown below. 
Considering the motion of the particle dX, and applying Newton's second 
law that a particle acted upon by an unbalanced force system has an 
acceleration in line with and directly proportional to the resultant 


of the force system, 


wae 
=z g =) ax = A dX = simplifying 
at 
oy (1) 
ce OSD. 2 (7 2D) 
oS Oe 3 ee 
and if E is constant, 
aD. EF OD 
at? PP axe 
Introducing a constant c of a valve such that of = —" this reduces 


to the familiar form of the wave equation, ref. (1h): 


ome c2 35D 
See Sx? 


Thus, it is established that equation (1) is a form of the wave equa- 
tion. It is desired to modify it to account for boundary conditions 


expected in a pile. This can be done in the general case as follows: 


9p Smee t ® , (2) 


with R representing resistance to driving. Equation (2) can be con- 
verted into a difference equation for numerical solution by standard 


methods (ref. 2h): 
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where E = and D(x,t) is the displacement of point x during time 
interval t. 

It is convenient at this point to introduce a simplified 
terminology which will facilitate comparison with model equations to 
be developed later, and to correlate with existing literature on this 


subject (ref. 28). 


Let D(x, t+ At) = D, 
Pex, 2), on 
mem) * 4, 

K 4 AX = = Ki 
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If the acceleration of gravity, g, is expressed in feet per second, 


then 


re 2 
Dy = 2dq - a + TBE Kaya = dy Mea = (Oy = dpe Ky = Ba (3) 


For purposes of numerical calculation, it is more convenient to use a 
mathematical model consisting of weights and springs such as is shown 
in figure (2). The following notation will be used and is consistent 
with that used in equation.(3) above and in reference (28). 

1, 2, 3, ...m-1l, m, .. . are subscripts designating par- 
ticular weights and associated springs, 

a, 2, 3, «. -. n-l, n, ... are subscripts desipnatime 
specific time intervals. Zero designates the initial instant. 


The following apply to any time interval n: 


D, = displacement of weight m measured from its initial 


position in inches 


Cn = compression of spring m in inches 

Fn = force exerted by spring m, in pounds 

Z_, = net force acting on weight m, in pounds 
NO = velocity of weight m, in feet per second 


The following apply to the preceding time interval n-l: 


d., = displacement of weight m measured from its initial 
position in inches 


c= compression of spring m, in inches 


m 

f{, = force exerted by spring m, in pounds 
ingles velocity of weight m, in feet per second 
z, = net force acting on weight m, in pounds 
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The following are usually constants: 
Wa ™ magnitude of weight m, in pounds 


= external force or resistance acting on weight mn, 


Ka = spring constant for spring m, in pounds per inch 
Rn 
in pounds 


= acceleration due to gravity (32.2 feet/second) 


Young's modulus of elasticity in pounds/square inch 


> tH 09 
n 


= cross sectional area of pile in square inches 


os 


= unit length in inches 


4’ t = time interval used for numerical calculations, 
in seconds 


rT, = critical time interval for spring m, in seconds 
TT. = mini j ec 
min Minimum value of T,, in seconds 


= over-all "factor of safety" = [oat 


Tin 


g 
a = individual factor of safety for springm = T,/At 


The numerical calculations will be a step-by-step process in 
which the five variables, D,, ©, FL, Z and V,, will be calculated 
for each weignt or spring in each successive time interval. Formulas 
for these five variables are developed by Smith (ref. 28) as follows: 

Dee, + yeleAt) (4) 

The coefficient 12 is included since Be and qd. are expressed in inches 
and Vj, in feet per second. 

Considering C, next, and referring to figure (3), let the dashed 
Squares represent the initial positions of weights m and m+ 1, and let 


the solid squares represent their positions at the end of time interval 


n, The initial length of spring mis &, and at the end of time interval : 
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it oul. D_ and Di.) are the displacements of weights m and m+ 1 


respectively, during time interval n, Then, 


Go, = £-£' 
but Le 7 = $' 7 De 
= Dn = Dae =L-£ 
therefore, Cc, = D-Day (5) 


ine ferce a is directly proportional to Cr and to the spring 
constant K., so that 


F, = OK (6) 


From figure (2), weight mis acted upon by springs m - 1 and 

m and by the external force Rs . the accelerating force is given by 
om = Fi 7 By 7 By (7) 

From Newton's second law that change in velocity is equal to 


force times time divided by mass, 


7 ate 
v7 %,* 2Z_( i ) (8) 


It remains to be shown that the model of figure (2) and equa- 
tions (i) through (8) are equivalent to the wave equation. This is 
done by using equations (1) through (8) in the following manner: 


Rewriting equation (8) for the previous time interval gives 


e+ 7) (8a) 


Similarly, equation (7) for the previous time interval is 


Zm = ie = fin = BR, (7a) 
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Substituting equation (7a) into equation ~ gives 


ie 5 + fe = £,- Rp ) Ate We 699) 


Substituting equation (9) into equation () gives 


De * by * [08 * Cg ~ f= Be) S88] 2208 


2 
12 at 
or D, = dy + 12 atv™ + ae [fn-1 - fm - Pm (ha) 


But equation (1) for the previous time interval is 


d, qe + v (AOR) 


=a 
ox wee (4b) 


Substituting (lb) into (ha) yields 


2 
= 4d, + le at [ny ‘i =e StE or -f,- Rin | 
m 


Substituting equations (5.) and (6) for the previous time interval 
into the above equation and a abe gives 

Dj, = 2dy - ot + AE © (due ~ dm)Kmea - (dm - yea )Km - Fon] 
which is equation (3). Thus equations () through (8) are equivalent 
to the wave equation, and the mathematical model of figure (2) has 
bee shown to be a valid one. 

Of course, solving the partial differential equation by a 
numerical integration process is an approximate one, but the results 
should be accurate within five per cent (ref. 28) provided that the 
time interval chosen is small enough and is coordinated with the length 
of pile selected to be represented as a block and spring in the mathe- 


matical model so that a stable set of calculations is produced. This 
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matter will be discussed in more detail later, but experience indicates 
that 5-to 10-foot sections of a pile may be so represented with a time 
interval of from 0,00025 to 0,00033 seconds, depending upon the type 

of pile material. As to the numerical method, it may be considered 
sufficiently accurate in light of present knowledge upon which the 


many factors influencing any pile driving operation are evaluated. 


B, Mathematical Model Representation 
That a pile may be treated mathematically as a series of 


weights and springs not only facilitates the computations but also 
assists in visualizing the problem, For that reason, the remainder of 
this paper will deal principally in terms of the mathematical model of 
a pile. It is, therefore, important that this method of pile repre- 
sentation be thoroughly understood at the outset. The following sec- 
tions will adhere closely to reference (28). 

The Hammer-Ram is usually a short, heavy, rigid object which 
may be represented by an individual weight without elasticity. In 
figure (4) the first weight W, represents the ram. 

The velocity of the ram at the instant of impact is needed to 
start the numerical calculations. Ordinarily the rated foot pounds of 
energy of the hammer is given by the manufacturer. The efficiency is 
sometimes given, and sometimes must be assumed. From these data the 
velocity may be calculated at impact by the following formula: 


velocity at impact __ mated ener ft-lb) x efficiency x 64. (10) 
(ft/sec ) weight (1b 
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The Capblock is a short, springy object of wood, plastic, or 
similar material which is comparatively light, and which may there- 
fore be represented by a spring. In figure () this spring is 
identified as K,. The form of the stress-strain diagram (or the 
hysteresis loop) that is produced as the capblock is suddenly com- 
pressed and then allowed to re-expand is assumed to be as shown by 
figure (5). Compression occurs along the line AB whose slope is 
determined by the elastic constant Ky of the capblock Restitution 
occurs first along the line BD, and then, because the capblock cannot 
transmit tension, it is completed along line DA. ABDA is the hyster- 
esis loop. The computations are made so that 

Area BCD _ fe 2 2 Energy output 
Area ABC il: nergy input 
where e) is the coefficient or restitution of the capblock. This 
relationship may be derived in the following manner: 
Let Pam be the weight of the ram 


W be the weight of the capblock 


cap 
Vram be the velocity of the ram before impact 
Veap be the velocity of the cap before impact 


v! be the velocity of the ram after impact 


ram 
vee be the velocity of the cap after impact 
Then, from Newton's laws, 


aie ca 


W V. ia ram ‘ram cap Yeap 


ram Yram ~ “cap Yeap 


and by definition, 
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i — 
If Veap = Yeap = 0, 
' 
Vv 
then e, = ram 
V 
ram 


2 
energy output (1/2) (Wran/8)(Veam) . (€1 Vpam)* 9 


rey input ~ al 
a = Ya) (weee/e) Gee) pepe 
From figure (5), 
area BCD _ ,2 . energy output (12) 


area ABC i energy input 


During compression the spring force may be computed from equation 


(6), 
F, = 0K (12) 


Inelasticity of the capblock must be considered during restitution, 
however, From equation (11) and figure (5), 


et = (1/2) BC x CD _ CD 


1 (1/2) Bx iC 
but, er) Co 
and CD = ef C 
1 1 max 


by similar triangles ADE and BCD, 


AE _ AD 
BC CD 
C - CD 
7 a elt BC 
CD 
aes a aia 
and AE = BC (—5—) 
e 
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Writing the equation for line DB gives 








P BC Baa j-e2 
1 7S % - B (4) 
CD eT 
and substituting for CD, 
aee 2 
e BC came Tol 
FS none Cy ~- BC ( 2 ) 
1 1 max iL 
Since BC « Cy aa KS : 
2 
FP Cy max KC, _¢ K Gat) 
9 1 max “1*",2 
e- CG L 
1 1 max 
which reduces to 
KC 
Ly 1 
a 1 


Calculations involving the capblock will, therefore, utilize 
alternately equation (12) during compression and equation (13) during 


restitution. © should not be thought of as a constant, but as 


1 max 
the maximum compression occurring in spring Ky before restitution be- 
gins. It often happens that recompression occurs and a new value of 
Cy max LS reached before restitution begins again. This is illus- 
trated by figure (5A) where the stress-strain diagram begins at A, 
goes to B during compression following equation (12), goes to C 
during restitution as governed by equation (13), at C recompression 
begins, goes to B following equation (13), continues to D by equa- 


tion (12), and thence to E by equation (13). As previously described 


the return to A is along EA because no tension can be transmitted. 
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That recompression can occur may not be apparent readily, but a 
feeling for it may be acquired in solving a problem or two manually, 
As will be noted later, equations (12) and (13) are used in 
the case of a cushion block, but with the subscripts changed to 2. 
Few tests have been made to determine the elastic character- 
istics of capblocks under impact conditions. Raymond International, 
Inc., in conducting a small number of tests found that the character- 
istics of a wooden capblock vary during driving, and concluded that 
in order to be on the conservative side in computing pile penetration 
per blow, a hardwood capblock with the grain vertical and of an 
original six-inch height, may be assumed to have the following 


characteristics: 
Spring constant K = 20,000 A (ib/inch of compression) 


where A is the horizontal cross sectional area in square inches. The 
coefficient of restitution was 0.5. The tests also showed that a 
Micarta capblock (Nema Grade "C") of 12 inches' height had a spring 


constant given by the following relation: 
K = 5,000 A 


and the coefficient of restitution was 0.8. 

Pile cap, or Follower, or Helmet--like the ram of the hammer, 
the pile cap is ordinarily a short, heavy, rigid object that can be 
represented by a single weight without elasticity, such as Wo in 


figure (4). If the pile cap is long and slender, as is the case when 
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used as a follower to drive piles below water or the ground, it is 
better to represent it by a number of weights and springs as shown 
in figure (6e). In this case elastic constants must be computed by 
the formula 


Ke = (Us) 


where A is the cross sectional area in square inches, E is the 
modulus of elasticity in pounds per square inch, and 4 is the unit 
length in inches represented by one spring. 

Cushion Blocks--(called "head packing in the U.K.)--in figure 
(4) springs K, to Ky inclusive represent the elasticity of the pile 
itself. However, if a precast concrete pile is being driven, a cushion 
block must be used under the pile cap W, to protect the concrete from 
shattering. In this case, spring K, would represent the cushion in 
combination with the first spring of the pile. Figure (6c) applies 
to the cushion block as well as the capblock, and they are treated 
mathematically in the same manner. Dynamic tests similar to those 
described for capblocks show that heinch-thick pine boards with the 
grain horizontal as used on top of a precast pile to distribute the 
blow evenly may be assumed to have a spring constant equal to 3,480 A 
and a coefficient of restitution of 0.5. Of course, the character- 
istics of the wood vary’ during driving, but use of the above 
relationship is on the conservative side when used for calculating 
penetration of the pile hammer blow. 

The Pile will be either of wood, concrete, or steel. All of 


these materials have traditionally been treated structurally as behaving 
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elastically within the working load ranges. Because the pile is 
compressible due to its length, it is subject to wave action under the 
blow of a hammer. This wave action may be analyzed by dividing the 
pile into short unit lengths of 5 to 10 feet. The weight of each 
length is noted as WF through Wijo, and the elasticity of each length 
ky an individual spring as shown by Ko through Ky in figure (ip ae 
the pile is of uniform section, the weights and springs representing 

it are identical. If the pile is tapered or is a composite pile, or 

1s in any way of irregular cross section, then the weights and springs 
are modified to represent the pile under investigation. That com- 
paratively large unit lengths, such as 10 feet, give desired accuracy 
in computing the action of impact waves may not be clear, Although 
water waves are transverse waves and the waves in a pile are longitudinal 
waves, water waves nevertheless can be used as an analogy to illustrate 
the principles involved. 

If a long strip of flexible material were allowed to float on 
the surface of a body of water, it would follow the wave action. If, 
for purposes of mathematical analysis, this flexible strip were rep- 
resented by a number of short floats of rigid construction but con- 
nected by flexible links, the mathematical model would appear as in 
figure (7). This representation would be an approximation, but 1t 
would involve negligible error because the small floats would ride 
the waves almost exactly like the flexible strip. If, however, the 
rigid floats are made comparatively long and they approach or exceed 
the length of the water wave, then the model would be as shown by 


figure (8). Obviously these long floats cannot follow the wave form 


closely. 
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Applying this analogy to a pile beirg hit by a ram, the 
importance of dividing the pile into unit lengths considerably shorter 
than the stress wave produced by a hammer blow can be appreciated. 
Fortunately, a pile driving impact usually produces a fairly long 
wave form, and unit lengths of the order of 5 to 10 feet will produce 
acceptable accuracy. In the special case where the exact form of 
the impact wave is being investigated, a smaller unit length of 1 or 
@ feet may be advisable. Under these conditions, it may also be ad- 
visable to divide the ram into a series of weights and springs as is 
done with the pile. For normal pile driving analysis, such refinements 
are not considered justified. 

There is a very important relationship between the unit length 
of pile selected and the time interval chosen for the calculations. 
the time interval to be used in the calculations is a function of the 
weight of the unit length selected and of the elasticity and density 
of the pile material. This relationship is derived and presented 
later, and governs the upper bound of the time interval to be used. 

On the other hand, choice of too small a time interval results in 
excessive computation for the accuracy required. 

The spring constants are determined by the formula given for 
K under "Capblock" above. 

The Ground--as indicated by figure (4) provision. is made for 
ground resistance to be applied on each section of the pile. In 
order to make a calculation using the wave equation, an assumption 
must be made as to the amount of resistance which will be offered by 


the ground at each of the pile sections. Furthermore, a stress-strain 
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relationship of the soil and viscous damping must be considered. For 
the purpose of this section, it will be assumed that these factors 

can be evaluated so that the method of their inclusion in the analysis 
may be developed. 

Resistance at the pile point--according to Chellis (ref. 5) 
the ground compresses elastically for a certain distance and then 
fails plastically with constant ultimate resistance Rie This concept 
is illustrated by figure (9). The elastic deformation will be called 
Nquake," and is denoted by the letter'Q.! The plastic deformation will 
be referred to as "permanent set," and is identified as 's's in figure 
(9) it is the distance AB or OC, There is a factor of importance 
which has not yet been included, and this is the rate of stress 
application. It is clear that the ground will offer more instanta- 
neous resistance to a rapid advancement of the pile tip than to a slow 
advancement. Viscous damping, therefore, will be included in the 
analysis to provide for the effect of rate of penetration. 

As has been shown, the wave equation calculation will give the 
instantaneous velocity of the point of the pile in any time interval. 
If J is designated as a damping constant and Ve the velocity of 
the pile point, then their product, Jv,, can be used to increase or 
decrease ground resistance at the point of the pile so as to produce 
damping. The instantaneous damping resistance at point x, for example, 
in figure (9) is Jv sR This value would be used for Ryo for a par- 
ticular fms interval applied as shown by figure (4). This damping 


resistance is a resistance to driving the pile and does not, of course, 
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contribute to the bearing capacity of the pile in supporting static 
loads. 

Resistance Along the Sides of the Pile--side resistance is 
calculated in the same manner as at the pile point, except that a 
different damping constant J’ is used, and it applies to Re through 
R in figure (4). As the pile is driven, the soil under the pile 
point is displaced and caused to flow aside very rapidly. On the 
other hand, the soil at the sides of the pile, as it is being dirven, 
is not displaced nearly as much, It may be assumed, then, that J! 
should be less than J. 

Distribution of Ground Resistance Between the Pile Sides and 
Point--the wave equation computation will utilize any distribution 
chosen. The problem of what distribution to choose will be discussed 
at length later, 

Figure (10) shows the model representation of the ground 
resistance alorg the sides of the pile. It must be remembered that 
ground quake is an elastic, or recoverable, deformation and that the 
set associated with it is of a temporary nature. The deflection of 
the spring Kk must, therefore, be limited to a maximum value less than 
or equal to Q, the ground quake. As the maximum allowed deflection 
of spring K is reached, provision is made for plastic deformation 
of the ground, es in the normal downward direction, or in an upward 
direction in the case of pile rebound. Figure (11) illustrates the 
way in which D) is handled in the calculations. It may be noted that 


De does not take on a value until D, exceeds Q. From that point on, 


Dt cannot be less than D, - Q or more than D, + Q. 
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The value of the ground resistance on a particular pile sec- 


tion is computed by the equation 
R, = (D, = DI) KI (1 + d'vq) (15) 


The point of the pile is dealt with in a similar, but slightly 
different, manner. Only downward plastic ground movement Dy is con- 
sidered. The maximum value of Dy is the permanent set of the pile 
which is designated as s . Figure (12) illustrates the way in 


which the calculations are handled. The value of the ground re- 


Sistance at the point is computed by the equation 


R, = (D, - DS) K, (1 + Jy,) (16) 


C, Illustrative Problem 

In order to bring the theory more sharply into focus, a 
hypothetical problem will be set up and solved manually using the 
general format described by E. A. L. Smith (ref. 27). 

Let figure (13) represent a steel pile 30 feet long. The 
weight of the ram is 5,000 pounds, and its velocity at impact as 
determined by equation (10) is found to be 10 feet per second. The 
capblock is hardwood, has a coefficient of restitution of 0.5, and 
its spring constant as determined by equation(5) is 2,000,000 pounds 
per inch. The pile cap weighs 500 pounds, ad is not able to transmit 
tension. No cushion block or head packing will be used, so K, will 
be the elasticity of the first pile section whose coefficient of 
restitution will be taken as 1.00. The pile will be divided into 


10-foot unit lengths, each of which weighs 500 pounds. The spring 
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constants K,, K3, and K), for the pile as computed by equation (6) 
are each taken to be 4,000,000 pounds per inch. Ground quake will be 
0.10 inch. The ultimate ground resistance has been determined to be 
220,000 pounds. It is distributed as follows: 10,000 pounds on each 
of the two top sections of the pile and 200,000 pounds on the bottom 
10-foot section of pile, Dividing the ultimate ground resistance by 
ground quake will give the soil spring constants as follows: 

! 


K, = i = 100,000 pounds per inch 


and Ke = 2,000,000 pounds per inch 


The time interval for the calculations will be 1/000 seconds. The 
damping constants J and J' will both be taken as equal to 0.1 for 
purposes of this example, although,as noted before, J! may be smaller 
than J. 

The calculations may then be performed by slide rule and 
recorded in tabular form as shown on the following pages, Although 
Space does not permit here, it is recommended that when making such a 
table a single large sheet be used. In referring to Table A, on each 
line values are first computed for D, ; Ds; D3, Dy Des then C,, Cys 
C3, ©, Cos etc., until all computations have been completed. The 
next line is approached in the same order, and so on, 

It may be of ‘interest to note from Table A that Cj reaches a 
first maximum at time interval 8; restitution takes place from time 
interval 9 through 22. Recompression occurs from time interval 23 


through 27 at which time a second C is reached. Restitution 


1 max 


continues until the end of the calculations, 
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C, = 
Di - Dy 


0 
0,03 
0.057 
0.078) 
0,093, 
0.1025 
0.1073 
0.1092 
0.1096 
0.1091 
0.1081 
0.1069 
0.1058 
0,1053 
0.1051 
0.109 
0.103 
0.1033) 
0.1022) 
0.1016), 
0.10238 
0.10479 
0.10855 
0,11266 
0.11687 
0.12095 
0.12h3 
0,12612 
0,12580 
0,1229h 
0.11925 
0.1168) 
0.11697 
0.11913 
0,121h2 
0.12185 
0.11971 
0.11586 
0.1127 
0.11129 


TABLE A 


Ce. 


(2x10°)(C, ), 


except during 
restitution 


60, 000 
11h, 000 
156, 800 
186, 800 
205, 000 
21h, 600 
218, L,00 
219, 200 
215, 200 
207,200 
197,600 
188, 800 
18h, 800 
183, 200 


181, 600 | 


176, 800 
169, 120 
160, 320 
Sco 
161, 440 
180, 720 
210, 800 
225. 320 
233,740 
ahi. 900 
28, 600 
252, 2h0 
2h9, 680 
226, 800 
197, ” 280 
178, 000 
179, O40 
196, 320 
21, 6,0 
218, O80 
200, 960 

170, 160 
113, 040 
133,600 


0 
- 60,000 
-11);, 000 
-156, 800 
~186 , 800 
-205, 000 
-21h,,600 
-218 , 1,00 
-219, 200 
=215, 200 
-207, 200 
-197, "600 
~188, 800 
-181, 800 
a1163) 200 
~181, 600 
-176, 800 
~169, 120 
~160, 320 
Co. 520 
~161, 0 
S1CO. 720 
-210,800 
-225, 320 
~233, * 7,0 
~24)1, 900 
~2):8, 600 
~252. 20 
~2):9, 680 
~226, 800 
eros 280 
~178, 000 
-179, 040 
~196, 320 
-21, 6,0 
~218, 080 
~200, 960 
fo 160 
Eis oho 
~133,600 
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0, 32506 
0.33792 
0. 34833 
0. 35652 
0.36334 
0.36897 
0.37361 
0.37781 
0. 3823h 
0.38796 
0.39475 
0.40127 
0.40556 
0.4065 
0.40) 
0.40136 
0.39910 
0.39835 
0.39835 
Ono fey 
0.39288 


TABLE A (cont'd. ) 








12, 000 
ho, 800 

8), 920 
13h, 280 
178, 000 
208, 320 
222. 800 
226, 200 
Paar "800 
220, 680 
216, 680 
210, 1480 
202,120 
193, 000 
186, 000 
183, 600 
187, 200 
197,920 
212,120 
226,680 
239,200 
29, 760 
25h, 360 
251, 000 
1,1. 8,0 
229,60 
205. 400 
232,320 
23,680 
2,8, 320 
239, "080 
218, 61:0 
197, le 
186, 8h0 
185,440 
196,680 
203,200 
190, "760 
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60, 000 
102, 000 
116,000 
101, 880 

70, 720 

36, 600 

1? "080 

-3,600 
-11,000 
~16, 600 
-23,080 
-27,880 
Soe 680 
-18. 5920 
-11, *}:00 
- 9. 200 
ere 4.80 
~26, 880 
-);2, 1,00 
aco 680 
-15,960 
~28 400 
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=20) 620 
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tee. * 600 
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It may also be noted that permanent ground displacement was 
not recorded until it exceeded the value assigned ground quake which 
was 0,10 in this problem. The time intervals at which these occurred 
were 

for D, , time interval 10 
Dy, , time interval 13 


De , time interval 15 


The calculations were terminated when all the velocities became nega- 
tive; effects after this are considered to be secondary and of little 
importance for the purposes of this study. 

A graphical presentation of the results is shown in figures 
(1h) through (16) inclusive. 

The reason that a time interval of 0.00025 seconds was chosen 
may not be clear, Generally, it may be said that the best time inter- 
val to use is one that produces a completely stable calculation. This 
means that when the results are plotted there should be no sharp 
Wiggles or reaks, If such irregularities are present, either an error 
has been made in the calculations or else the time interval chosen is 
too large. Once the unit lengths of the pile have been chosen, the time 
interval must be selected to correspond. Choosing too small a time 
interval would result in an excessive amount of calculation with only 
a very limited increase in accuracy. 

Each spring in the figure (13) has a critical time interval 
which is the time that is required for a stress wave to traverse the 


spring and its associated weight. Stress waves travel in both directions 
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at a speed of -VE/P where f is the mass per unit of volume. The 
following equations may be derived for the critical time interval for 
spring Ra which is called T,. 


For wave motion downward: 








a “Veith - lee 


For wave motion upward: 





yee yt 
m “Vioek- ~ 19.648 Kn 
The smallest value of T,, obtained by using these formulas is the 
critical one. About one half of the critical value should be selected 
as the time interval to be used in any particular problem so as to 
prevent instability from arising due to other factors not included in 
the above formulas such as quake, damping, and coefficients of restitu- 
tion. When dividing the pile into ten-foot sections for purposes of 
mathematical representation, the recommendations of Smith (ref. 28) 
have been found to be generally satisfactory: 

For steel piles -- 0.00025 seconds 

For wood piles == 0,00025 seconds 

For concrete piles -- 0.00033 seconds 
During some of the test cases run, instability resulted in using this 
time interval for wood piles, and recomputations were made using smaller 
time intervals. Instability can generally be detected by plotting 
the results as noted previously or by comparing the velocity of the 
pile cap and the pile tip with the velocity of the ram at impact. If 
either of these velocities exceeds twice the ram velocity at impact, in- 
Stability is probably present. The velocity check should not be relied 


upon alone, however, to establish the stability of the calculations. 
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IV, COMPUTER ADAPTATION 


A. General 

As shown by the illustrative example, pile driving may be 
analyzed by the wave equation uSing a numerical solution method and a 
Slide rule, ‘This approach is a valid one and will produce satis- 
factory results. The main problem, of course, is the amount of time 
and work required to make the calculations. The example problem was 
selected as being about as simple as possible, with convenient figures 
to deal with numerically, and still illustrate all aspects of the cal- 
culations. Nevertheless, it required about twelve hours of manual 
computation to get the final set. It should also be noted that an 
error made at the beginning of a solution will carry throughout the 
calculations. The usual problem encountered in practice would probably 
take several mandays to complete. 

Fortunately, there is a quicker and better way of making these 
calculations. Use may be made of high-speed digital computers which 
can solve up to 800 such practical problems in about 5 minutes, de- 
pending upcen the computer used, As a consequence, use of the wave 
equation for pile driving analysis has become a practical matter. 

Although it is not the purpose here to go into the details of 
programming, it may be well to mention some of the broader aspects for 
those totally unfamiliar with computers. For a computer to be effective 
in solving engineering-type problems, the problem-solving procedures 
mus t be presented to the computer in a language which it can "understand." 


The computer's basic language consists of elementary instructions such 
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as add or subtract. The problem-solving procedure must be translated, 
then, into instructions which the computer can obey. This translation 
is called programming and can be carried out entirely by a person, or 

the computer can assist in the process by use of a compiler. 

A compiler is a large set of computer instructions which can 
accept a problem-solving procedure, written in a form resembling the 
procedure, and produce from it the proper elementary machine instruc- 
tions that will solve the problem. The algebraic compiler system used 
for the programs here is called Fortran, which stands for formula 
translation. A procedure to be followed in solving a problem with 
Fortran is specified by a series of statements, of which there are 
several types. One type specifies arithmetic operations: another calls 
for reading a data card or printing results; a third specifies the 
sequence in which the statements are to be executed; and a fourth type 
is made up of statements which provide information about the procedure 
without themselves causing any action. All of these statements taken 
together form a source program. The programs included in the appendices 
of this writing are source programs. When punched on cards and put 
through the machine, it is translated by the Fortran system into simple 
sets of machine instructions which is called an object program. For 
those who wish to pursue the subject further, reference (21) provides 
in the space of 83 pages all the information necessary for a person to 


learn to develop programs of his own. 


B, Description of Programs Developed 


1. Basic Programs 
A basic program in Fortran language for the IBM 7090 computer 


was prepared, and it follows the method of calculation illustrated in 
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the example problem. This program requires that input data be supplied 
to indicate the number and size of the weights in the mathematical 
model, the time increment to be used, the spring constants of the pile, 
the spring coefficients of the soil, the coefficients of restitution 

of the capblock and cushion block (if used), the velocity of the ram at 
the instant of impact, the value of ground quake, values of point and 
side damping, whether or not the pile cap can transmit tension, and 
whether or not side resistance is present, Output of the basic program 
includes the input data for purposes of identification and checking, 
and then for each time increment the force and resistance applied to 
each block of the model together with displacerent and velocity of each 
block at that particular instant. The permanent set of the tip of the 
pile is also recorded. The prograin is arranged to stop the computer 
after 300 time-intervals or before if 

(a) the velocity of the pile cap exceeds twice the velocity the 
ram had at the instant of impact; 

(b) the velocity of the pile tip exceeds twice the velocity the 
ram had at the instant of impact; 

(c) the velocities of all the blocks are simultaneously negative 
or equal to zero. Velocities toward the top of the pile are taken to 
be negative. 

The basic program can be used to study stresses, velocities, 
accelerations, and permanent set of the pile as they vary with time. 
The program may be run with one or more sets of data at a time. This 
program in its entirety, together with instructions on how to use it, 


1S included as appendix A, 
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Having developed this program, the question was how best to 
use it in trying to determine if correlations could be achieved with 
available field data on driving records, load tests, and borings. Of 
the variables that must be provided as input data, the greatest un- 
certainty was with the values to assign ground quake, point and side 
damping, and distribution of resistance along the sides of the pile. 

2. Vary RU Program. 

Correlation was approached in two different ways. The first 
- method was particularly suited to the very few cases in which careful 
and complete field tests were reported and analyzed comprehensively. 
In these cases, reasonably accurate estimates of the variables except 
ground quake and damping factors could be made. Another program was 
developed for this approach which was basically the original progran, 
but automatically for the one set of input data varies ground quake, 
ultimate ground resistance, and damping factors over a wide range of 
values. The output includes for each variation of basic information 
indicated the ultimate ground resistance, the blows per inch at final 
set, maximum tension and compression at the head, mid-length, and tip 
of the pile, the time interval the calculations stopped and why, and 
the original input data. ‘his program in its entirety is included 
with instructions on how to use it as appendix B. 

3. Researcher Program. 

The second method of approaching correlation was more of a 
statistical approach, and it was designed to try and glean as much 


information as possible from published field data which were not as 
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comprehensively reported as data used in the first method of approach. 
References (15) and (16) were used extensively in this attempt, Certain 
pile sections, lengths, and driving hammers were chosen because more 
reasonably good data were available for them than for other types, 
Another computer program was developed around the original program and 
designed to systematically produce in quantity data from which curves 
of ultimate ground resistance vs. number of blows per inch could be 
plotted for various values of damping constants and ground quake and 
side resistance, This program, together with instructions for its use, 
is included as appendix C, It is designed to be run with one or more 
sets of data. Its output includes for each of the above automatic 
variations of input, the ultimate ground resistance, values of soil 
spring coefficients, the permanent set, the blows per inch at the final 
set, maximum tension and compression at the head, mid-length, and tip 
of the pile, the time interval the calculations stopped and why, and 
the original input data, Thus about eighty complete problem solutions 
are obtained for each set of data run. Ordinarily three sets of data 
can be run at a time in this manner and take no more than five minutes! 
IBM 7090 computer time. Manual solution would require at least one man- 
year for just one set of data, Running time on a 709 computer would 
be 20 to 25 minutes, but only about a minute if the 709) is used. 

4, Validation of Computer Programs 

Before using any computer program, great care must be taken to 
prove that it is correct and that all of its features behave as they 
are intended to. Making a manual solution is the only way to be cer- 


tain that the program is correctly written. To prove the original, or 
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basic, program, the data of the illustrative example presented 
herein was used as input, and the solution compared with the manual 
solution. A tabulation of the results of the two solutions for 
displacements of the ram in each time interval are found as Table B. 
It can be readily seen that the results compare as closely as one 
would expect considering that a slide rule was used in the manual 
computation. As a further check of the program, the illustrative 
problem presented by E. A. L. Smith (ref. 28) was computed using the 
basic program, and the results plotted. They compared in every 
detail with Smith's solution. 

As additional programs were developed around this basic 
program, they were tested with both of these illustrative examples 
and verified to insure that subsequent work would be based on sound 
programs, 

5. Computer Program for the Hiley Dynamic Formula 

The Hiley type formula is thought by many to be one of the 
better dynamic formulas in general use today. Although not accepted 
by all, it is becoming more popular in the United States and has been 
included in some building codes. In order to facilitate comparison 
of results obtained by the wave equation and those obtained by using 
the Hiley formula, a computer program was developed for it also. 
This program requires that input data be provided to give the hammer 
data; pile material; the pattern, either rectangular or triangular, 
of side resistance distribution; coefficient of restitution; weight 
of the pile; area of the pile at the top, mid-length, and bottom; 


the ground quake; the modulus of elasticity of the pile material; 
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Table B 


COMPARISON OF COMPUTER SOLUTION WITH MANUAL SOLUTION 
FOR ILLUSTRATIVE EXAMPLE 


Displacement of Block 1 


Time Interval Computer Results Manual Solution 


a 0,030 0,030 

2 0,060 0.060 

3 0.089 0.089 

hy © a7 0,346 

E 0.145 0.145 

6 Oean7a: 0.172 

t 0.197 0.197 

8 (0) zal O, 221 

9 0,25 0,25 
10 0,267 0,267 
ale 0,288 0,288 
12 O, 309 0.309 
13 0.328 O. 328 
Wy O.3h7 0.3h7 
15 0,364 0,365 
16 OABer O. 381 
17 0397 0.397 
18 0.412 Cai 
19 0,426 0.426 
20 O.))0 0.440 
on 0.453 0.453 
22 0.65 0.65 
23 0.476 0.476 
2h 0.85 0.485 
25 O.L94 0.95 
26 0,502 0.502 
27 0,508 0.508 
28 0.513 Oy 5a 
29 0.517 0,518 
30 0.520 0.520 
31 0,522 7522 
32 0,523 0.523 
Be 0.523 0.52h 
3h 0.522 0,523 
25 0,522 0,521 
36 0.517 0.518 
Bf C7583 O7515 
38 0.508 0.509 
39 0.502 0.504 


END OF CALCULATIONS 
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and the pile length in feet. The-program will automatically range 
over ultimate ground resistances from 20 to 380 tons in increments 
of 20 tons, writing out the set and blows per inch corresponding to 
each ground resistance as well as repeating the input data for each 
case. It computes for each value of the ground resistance the set 
associated with all the load carried by the point, one half of the 
load carried by the point, and none of the load carried by the point. 
Programming the Hiley formula may not be justified under usual cir- 
cumstances, but it was done in this case as a convenience. The 
program and instructions for its use are included as Appendix D. 

It should be noted that the Hiley formula is discriminatory to long, 
heavy piles and is subject to the same general limitations as any 


dynamic formula, 
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V. SOIL ENGINEERING CONSIDERATIONS 


As has been mentioned in the development of this method of 
pile driving analysis, values must be assigned to the dynamic soil 
properties of ground quake and damping. The distribution of forces 
acting on the pile tip and sides must also be evaluated to complete 
the information required for a solution. Furthermore, in order to 
approach correlation between dynamic driving resistance and static 
resistance, the manner in which frictional resistance along the 
sides of the pile acts in both the static and the dynamic cases must 
be explored. The subject of load testing must also be investigated 
insofar as the way in which it is done and the interpretation of 


results affect correlations, 


A, Side Friction Under Static Loading 
1. Granular Materials. 


Some recent work has been done on load transfer in the case 
of end bearinz piles in which two instrumented H piles were tested 
(ref, 12). ‘the piles were 1)BP89 and 1)BP117 and were driven ap- 
proximately )O feet through granular soil varying from sandy silt 
to sand and gravel, The piles were end bearing on weathered silty 
shale. In spite of the fact that the piles were short and rela- 
tively stiff, they transferred as much as one third of their load 
to the surrounding soil by friction. Figure 17 shows the soil con- 
ditions and driving records, 

Upward shear along the pile develops as the pile shortens 


elastically under the applied static load at the top of the pile. 
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This movement of the pile causes distortion of the surrounding soil 
resulting in development of a frictional force resisting the dis- 
placement of the pile. Generally, it would be assumed that the 
greatest side resistance would be developed near the top of the pile 
and would diminish with depth. This action may be visualized by 
thinking of the pile as a spring embedded in the ground with a load 
applied at the top of the spring. One might also expect a nearly 
linear relationship between load applied and load transferred as long 
as no relative movement occurs between the soil and the pile. How. 
ever, a loading is reached at which the soil begins to slip along the 
pile. This zone of shear failure begins at the top of the pile and 
progresses downward. If the pile tip never moved, the entire applied 
load would be carried by friction along the pile sides, but even in 
the case of progressive domward failure, the failure zone would 
never extend down to the tip. In the actual case the pile tip does 
move downward as the applied load is increased so that it accepts a 
portion of it. 

The test results reported in reference 1 were presented from 
a Viewpoint of over-all transfer of load to the soil by friction. 
Taking the data presented, and dividing each of the test piles into 
sections about 5 feet long, Table C has been prepared in an effort 
to show the variation in friction force on each section of the pile 
as the test load was increased, In each of the test piles, section 
1 is the top section. A study of this table will show the validity 
of the description of the interaction between pile and soil, i.e., 


the frictional forces have a tendency to develop more at the top of 
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Table C 


| FRICTION FORCE ON SECTION 
(IN TONS) 


Increment of Total 0-75 Tons 75-150 Tons 150-225 Tons 225-300 Tons 


| Applied Load Change Total Change Total Change Total Change Total 
(at 75) (at 150) (at 225) (at 300) 


Section of 1)BP89 
Pile and Soil 





Description 
1e-Sand & Gravel 8 8 2 10 5 15 hy 19 
2—-Sand & Gravel i if 2 9 hy 13 2 15 
3--Fine to Med,Sand 6 6 1 7 -2 5 3 8 
le-Fine to Med,Sand 0 0 3 3 hy 7 al 8 
S§--Sand & Gravel 5 5 2 7 =3 by 4 8 
6--Sand & Gravel hy h 0 by a 5 8 13 
7--Sand & Gravel 3 3 = 2 5 if 5 ae 
8--Sandy Silt, tip 
in weathered 
shale 1 iE 5 6 0 6 h 10 
TOTAL FRICTION FORCE 3h 48 62 93 
TOTAL FORCE ON TIP 1 102 163 207 
Section of 1)BP117 
Pile and Soil 
Description 
1--Sand & Gravel h hy <1 3 -3 0 0 0 
2--Sand & Gravel by hh i 5 -3 2 0 2 
3--Fine to Med.Sand 0 O yy ly -); O 2 2 
l--Fine to Med.,Sand 2 2 3 5 -5 O 1 i 
5e-Sand & Gravel 3 3 1 hy cS aly 11.5 28.5 
6--Sand & Gravel 3 3 1 hy 4 8 =e 6.5 
7--Sand & Gravel 6.5 6.5 3.5 10 -1 9 5 1h 
§--Sandy Silt, tip 
in weathered 
shale 0.5 0.5 1.5 2 hy 6 -l; 2 
TOTAL FRICTION FORCE 23 Bt 2 5h 
TOTAL FORCE ON TIP 52 113 183 2h6 
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the pile. It also shows that local failure extends over a wide range 

of total applied load, For example, as the applied load is increased 
from 75 to 150 tons, the friction forces on section 7 of the 1) BP 89 is 
decreased by one ton; from 150 to 225 tons' applied load,local failure 
occurs at sections 3 and 5. After this loading the granular structure 
of the soil seems to have been readjusted by the shearing forces and 
displacements so as to have an increased load-carrying capacity. 

In the case of the 1 BP 117 pile, as the applied load is in- 
creased from 75 to 150 tons, local failure has taken place on section 
l only. As the applied load is increased to 225 tons, failure continues 
on section 1 and is expanded to include sections 2, 3, 4, and 7. Upon 
further increase of the applied load to 300 tons, it is evident that 
failure has generally progressed down to and including section even 
though there is a very slight increase in friction force on the top 
half of the pile. Sections 6 and 8 are experiencing failure while 
sections 5 and 7 are attempting to assume the load. Even at this 
depth it is interesting to note that section 5, the uppermost of these 
two sections, has assumed an appreciably higher part of the load. 

Both of these piles were designed for 150-ton static loads. 
When test-loaded to 150 tons, it is noted that the general shape of 
side frictional resistance was triangular with the base of the triangle 
at the top of the pile and the apex toward the lower sections. In the 
case of the 1) BP 117 at a test load of 150 tons, however, the general 
shape of the frictional resistance acting along the length of the pile 
was generally of a rectangular shape down almost to the pile tip. Thee 


difference may be explained by considering the loads required for a 
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unit axial shortening (spring constant) of each of the pile shapes. 
From the expression K -= » and the fact that the length and modulus 
of elasticity of each pile is the same, their spring constants! ratio 
will vary as their area ratio or 3h./26.19. It can be seen that the 
1h BP 117 is 1.32 times as "stiff" as the 14 BP 89. Since the over-all 
dimensions of the piles, the applied load, and the surrounding soil 

are essentially the same, this relative "stiffness" must account for 
the difference found in the pattern of frictional resistances developed. 
As a consequence of this, it is interesting to note that the 14 BP 117 
was nearer total failure at a test load of 300 tons than was the 
lighter 1h BP 89 pile. 

In summary, the following factors may be important in con- 
sidering the static behavior of piles in granular material: 

a. The time which has elapsed since driving. In some cases a 
quick condition may have been caused during driving which will dis- 
appear as pore pressures dissipate, Relaxation may also occur. 

b. The general manner in which load is transferred to the 
soil may vary considerably and is generally dependent upon the "stiffness" 
of the pile in resisting axial deformations as related to the strength 
characteristics of the soil. The less "stiff" piles tend to give up 
their loads to the soil at the upper strata while the very "stiff" 
piles tend to carry their load by end bearing. 

c. Failure of piles in granular materials is often not too well 
defined in the case of friction piles. Small local failures along the 
pile may result in a decrease in load carrying capacity, but then the 


soil particles in undergoing shear may assume different packings and 
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become denser with the result that their strength has been restored 
or even increased over its original value. The pile load may be 
increased while these readjustments are taking place until a load is 
finally applied which is sufficient to push the pile through the soil 
continuously. With each of these readjustments the soil takes on 
different properties causing the pile to undergo variations in its 
bearing capacities. 

d. Unusual conditions which may be pertinent to a particular 


problem such as a quick condition or hydrostatic uplift. 


2. Cohesive Materials. 

As in the case of piles in granular material, some research 
work has also been done with instrumented piles in cohesive materials. 
H. B. Seed and L. C. Reese report in reference (2) that several 6~inch- 
diameter pipe piles from 20 to 22 feet long were driven approximately 
15 feet into a stratum of soft, saturated clay. Electric strain gages 
were installed on one of these piles. Test loadings to failure were 


made at various intervals after driving as follows: 


Interval After Driving Failure Load 
~ (pounds) — 


pounds 

3 hours 1070 
1 day 2800 
3 days 3200 
7 days 51,00 
1h days 5800 
23 days 6100 
33 days 6200 


This data show that the ultimate bearing capacity increased very 
rapidly with time for the first week, but that thereafter the increase 


was very gradual. 
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Most of the load was removed from the pile by side friction 

with only about 10-15% of the applied load reaching the pile tip. 

The increase in supporting capacity was reflected by an increase in 

the load transfer to soil at all depths as shown by Table D which was 
derived from published data (ref. 26). The pile was divided into 
three-foot sections with section 1 being at the top of the pile. From 
this table it can be seen that the resistance is generally increasing 
with depth for the first 9 to 12 feet of embedment, and then it tends 
to become constant with depth for the final 3 to 6 feet of embedment. 
It is interesting to contrast this behavior with that already described 
for granular materials. 

In this series of tests relationships were established for 
undisturbed soil. and fully remolded soil. Tests of soil samples taken 
next to the pipe wall showed that the strength loss caused by dis- 
turbance of driving the pile was about 60%to 70% of that which the soil 
would have lost due to complete remolding. The clay next to the pile 
wall subsequently gained strength, and at the end of the test period 
it had a strength of 60% higher than the strength of the undisturbed 
soil. This increase in soil strength appeared to parallel the decrease 
in excess hydrostatic pressure indicating a relationship between the 
two. For all values of time at all depths the failure of the friction 
pile did not occur in the soil, but at the interface between pile and 
soil. 

In another series of tests made by Vey (ref. ooo. with H piles, 
he found that in some cases the friction force reversed itself when the 


applied load was removed and held the stress in the pile. He also 
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Table D 


FRICTION FORCE ON SECTION 
(IN POUNDS ) 


Test 7--33 Days After Driving 





Increment of Total Q-2000 lbs. 2000-000 lbs. 4000-6000 lbs. 
Applied Load Change Total Change Total Change Total 
(for (for (for 
Section of 6" Pipe 2000) 4,000) 6000) 
Pile and Soil 
Description 
1--Soft sat. clay sae a0 0 110 50 160 
2--Soft sat. clay 290 290 270 560 150 710 
3--Soft sat. clay 480 480 1,00 880 h20 1300 
l--Soft sat. clay 1,80 4,80 520 1000 500 1500 
S--Soft sat. clay “75 175 525 1000 560 1560 
TOTAL FRICTION FORCE 1835 3550 5230 
TOTAL FORCE ON TIP 165 SO Te 
First Loading to O0-1070 lbs. 
Failure (3 hrs. Change Total 
after Driving) (for 
1070) 
le-sam soil O 0 
2-- 200 200 
3—— O O 
= 330 330 
5-- 280 280 
TOTAL FRICTION FORCE 810 


TOTAL FORCE ON TIP 260 
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found that adherence would increasé between the soil and the pile with 
time, and reconsolidation of the disturbed clay would increase its 
shear strength. 

The findings of these two series of tests regarding the effect 
of pile driving on clay cannot, of course, be assumed to be valid under 
all conditions. Tschebotarioff has shown in reference (3) that the 
remolding effect caused by driving pipe piles closed end in 100 feet 
of varved clay can be great. Remolded clay layers can lead to the 
development of negative friction with the result that end bearing 
piles may become overloaded, It is generally understood that displace- 
ment piles should not be driven in highly sensitive clays; therefore 
this part of the problem will not be discussed further. 

In summary, it may be said that the static behavior of piles 
driven in cohesive material may depend upon: 

a, The time which has elapsed since driving. In some cases 
the increase in bearing capacity for even short periods of time may 
be great; in other cases a loss of bearing capacity may result. 

b. The friction force acting on the pile tends to increase 
with depth. It is believed that the load transfer would decrease with 
depth if the shear strength of the clay were sufficient to accept the 
loads at the upper strata and if there were a motivating force avail- 
able to cause these upper strata to come into intimate contact with 
the pile. ‘Since neither of these conditions usually exists, the load 
transfer is found to generally increase with depth. 

c. Failure of piles in clay may be more dramatic than those in 


granular materials. In the case of a pile not resting on an extrenely 
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hard strata such as rock, once the shear strength or adhesion has been 
exceeded, the pile may begin to punch its way through the supporting 
strata and continue downward until the load is reduced or the pile 


encounters a firmer strata. 


B, Dynamic Factors Involved in the Pile-Soil Interaction 


The main obstacle in attempting to relate the penetration 
record obtained during driving of a pile to the static load carrying 
capacity is that of obtaining a reasonably accurate estimate of the 
dynamic resistance to penetration between the pile and soil. The prob- 
lem is the same whether using the wave equation or a dynamic pile 
driving formula. 

The dynamic resistance to penetration is related to one or 
more of the following factors: 

a. The static resistance. 

b. The increase in strength observed to occur in the soil 

due to rapid loading. 

c. Frictional resistance between the pile and soil. 

d. Viscous resistance developed as the pile penetrates the 

soil. 

The basic strength of the soil, whether due to cohesion or 
internal friction, constitutes the static load capacity of the pile 
and soil. At least, this represents the limiting value since the pile 
may not develop the full shear strength of the soil when acting as a 
friction pile. The value of the static resistance can be thought of 


as the base value of the dynamic resistance in the idealized case 
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where the strength properties of the soil are not changed by the 
process of pile driving. This is rarely the case, and usually the 
driving of the pile acts to decrease or increase the natural strength 
of the soil. 

Several investigators have studied the effect of rapid loading 
times on the compressive strength of the soil (ref. h, 11, 25, 37), 
and some discussion of the results of such investigations appears 
below. It seems logical to assume that in the early time intervals 
of pile penetration most of the dynamic resistance is a result of 
this type of phenomenon, the magnitude of which would be directly 
related to the static strength of the soil. 

Frictional resistance in nearly all cases develops between the 
pile sides and soil as a result of the penetration of the pile. The 
nature and magnitude of this frictional resistance depends upon the 
type of soil, moisture content, type of pile and pile material. If 
the adhesion or friction acting between pile and soil is greater than 
the shear strength of the soil, the frictional resistance may be de- 
veloped within the soil itself. The dynamic frictional resistance 
would be of about the same value as the static resistance since the 
dynamic coefficient of friction between the pile and soil materials 
is about equal to the static coefficient of friction (ref. 2h). Fric- 
tion as a dynamic factor will. be neglected since there is some doubt 
as to the extent that actual relative movement takes place between the 
pile and soil at the pile-soil interface. 

No published work has been found that deals with viscous 


forces as developed by penetration into soils. One investigator 
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attributes part of the influence of high rates of strain on the increase 
in compressive strength of soil to viscous action (ref. 37). Since 
viscosity is defined as a property of internal friction (ref. 2), it 
will be considered as included as one of the factors contributing to 

the increase of compressive strength experienced by a soil subjected 

to high rates of strain. 

During the driving of a pile, it is assumed that the resistance 
developed at the tip and along the sides can be expressed in terms of 
the resistance developed under static loading. This is to say that 
the dynamic resistance to driving is composed of the static resistance 
plus an increment of resistance that develops under dynamic loading 
and is expressed as a percentage of the static value. The problem is 
to determine what form the relationship between static and dynamic 
resistance takes and how to quantitize it. 

Another problem, which will be taken up later, is how much of 
the static resistance is present during driving. It cmnot be generally 
assumed that all of the static resistance is present during driving, 
and Chellis points out ways in which this may vary (ref. 6). On the 
other hand, dynamic resistance may develop in soil strata that offer 
only very small or negligible static resistance. 

In the computer solution, the total resistance to driving at 


the pile point is represented by the expression: 
R = (D, =D) kK’ (1 + Jv) (15) 
where CD... = iy is the "elastic" strain of the soil with 


a limiting value of Q, 


K' is the spring constant of the soil, 
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v is the pile tip velocity at any instant, and 
J is the damping constant, 

This relationship is shown by Smith (ref, 28), 

Two questions arise with respect to this mathematical representa- 
tion of the dynamic resistance to driving: 

(a) Does this equation approximate the actual form of the pile 
reSistance; i.e., does the resistance-time relationship approximate 
that of the actual case? 

(b) What values of J should be assigned for different soil and 
driving conditions? 
Smith proposed the use of a value of 0.15 for J; however, since the 
purpose of his published work did not include the correlation of re- 
sults with load tests, he gave no basis for the use of this value. Some 
preliminary attempts to correlate wave equation results with the results 
of pile load tests and driving records using values of 0.15 and 0.05 
for J and J! indicated that this value for J was not appropriate far 
clay soils. 

In an attempt to answer the questions posed above, a search of 
the literature was made for any experimental work that might lead to 
a confirmation of the nature of the dynamic resistance as well as toa 
determination of more appropriate values for the coefficients J and J!. 

Several research programs over the past 15 years have investi- 
gated the influence of time of loading, or strain rate, on the compres- 
Sive strength of soils. All such projects have used rapid loading in 
triaxial and unconfined compression tests, or dynamically loaded footing 


models as the method of investigation. The application of results fran 
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such tests to pile driving must be approached with caution and should 
be considered only as indicative of the general nature of the dynamic 
resistance to pile driving : and of the general order of magnitude for 
values of the coefficients of damping. 

Before describing the tests made by individual investigators, 
several terms will be defined: 

(1) Time of loading: The time from the beginning of a test 
to that at which the maximum force is achieved. Figure (21) shows 
this quantity graphically. 

(2) Strain rate: The rate, expressed in per cent per second, 
of strain imposed on the sample at the time the compressive strength 
is reached, 

(3) Transient load: A dynamically applied load in an uncon- 
fined compression or triaxial test that reaches a maximum in times of 
less than 30 seconds and as short as 0.005 second. 

(4) Static test: For comparative purposes, a test that is made 
with a time of loading of ten minutes or more, unless otherwise indicated, 

(5) Strength ratio: The ratio of the compressive strength for 
a given time of loading to the strength for a static test. 


(6) Strainerate effect: Same as strength ratio. 


2. Work at Harvard 
In a test program sponsored by the Panama Canal in 19)8, 
Casagrande and Shannon made a series of static and transient unconfined 
and triaxial tests on soft to medium clays and a soft rock (ref. )). 


Triaxial tests of the vacuum type were also made on a dry sand. Times 
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of loading for the transient tests ranged from 0.01 seconds to 30 
seconds. 

In all the tests on clays, the transient compressive strength 
was greater than the compressive strength as determined by static tests. 
The transient strengths, taken at a time of loading of 0,02 second, 
were from 1.4 to 2,6 times the static compressive strength. The tests 
also showed that the increase of the fast transient compressive 
strength over the static strength is greatest for specimens at the 
highest water content and least for the specimens at lowest water 
content. (Figure (17-A) is a summary of the average results from the 
clay tests. 

Vacuum-type triaxial tests were performed on a clean, medium 
sand with a void ratio of 0,62 and at a minor principal stress of 
2,3 kg/cm. The test results indicated that the compressive strength 
of a dry sand in fast transient tests is about 15% greater than the 
static compressive strength. 

3. Work of Whitman 

Whitman measured the increase in compressive strengths of 
solls using a triaxial machine, and testing at high rates of strain 
(ref. 37). Table E shows his results for several types of cohesive 
soils. The last column of this table (S-R) is the ratio of the 
failure load at a strain rate of 1000 per cent per second to the 
failure load at a strain rate of 0.03 per cent per second. Itis 
observed that the strain-rate effect varies from 1.3 to 2.0 for these 
soils. In his tests, Whitman achieved times of loading as short as 


.005 second. Confining pressures up to 100 p.s.i. were used, 
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Table E 


INCREASE OF COMPRESSIVE STRENGTH 
IN COHESIVE SOILS DUE TO HIGH STRAIN RATES 


(from ref. 37) 


DESCRIPTION Pol, Pel. ly le cueeeeeen 
1. Normally consolidated, sensitive 
ocean sediment, undisturbed 63 49 92. O53 azao 
2. Very plastic clay, remolded 27 38 h8 7 Sages 
3. Plastic clay, remolded 17 am 16 ieee? 
4, Medium soft, slightly sensitive 
clay, undisturbed 2h 26 a TO Meee 
5. Slightly plastic clay loan, 
remolded 23 22 ol 3 ee 
6, Plastic clay, remolded 27 38 aise ey 
1. Moderately sensitive silty clay, 
undisturbed ee 22 35 22) Gao 
8. Impervious compacted fill a7 bal ale 25 1.8 
9, Tough compacted fill yi 21 26 35 2.0 
10. Stiff, dry clay, undisturbed 23 30 20... 250 Win 


in comparing results of tests on clay with and without confining 
pressure, Whitman observed two different time effects. The first effect 
can be likened to that encountered in an extremely viscous fluid, and 
is associated with a continuous, plastic deformation of the soil. The 
second effect is associated with the formation of discontinuities such 
as shear planes and splits, md occurs when the soil sample is tested to 
failure without (or small valves of) confining pressure. From these 


results, he concluded that when cohesive samples were subjected to a 
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lateral pressure, the strain-rate effect was the result of viscous 
phenomena, and was the same regardless of the magnitude of strain. 
The same effect occurred in very plastic soils regardless of confining 
pressure. 

In dry or moist sands an increase in compressive strength of 
10-15 per cent was observed. Only limited tests were conducted with 
saturated sands, and figure (18) shows results for saturated Ottawa 
sand and a well graded sand. Whitman observed that the strain~rate 
effect in saturated sand must result from differences between the pore 
water pressures in slow and rapid tests. Whensand is tested at 
constant volume, pore water migrates to the central part of the sample. 
As speeds of deformation increase, more energy must be expended to 
overcome resistance to the flow of water. 

4, Work of Seed and Lundgren 

Seed and Lundgren investigated the effect of transient loadings 
on the strength and deformation characteristics of saturated sands 
(ref. 25). Their tests were also made with the triaxial apparatus, and 
three different categories of tests were conducted: 


a. Static tests in which the loading times were about 
10 minutes; 


b. Slow transient tests with loading times of about 
lh seconds; 


c. Rapid transient tests with loading times of about 
.O2 second. 


In the rapid transient tests a constant rate of deformation of 0 inches 
per second was used, All tests were conducted with a confining pressure 


of 2 kg/sq.cm., and the samples were tested at various void ratios. 
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Results of the tests on dense, fine sand indicated that the 
compressive strength of a sample subjected to rapid transient loading 
was about lO per cent greater than that obtained by static drained 
tests. Due to the high rate of loading the rapid tests were actually 
undrained tests, since drainage could not occur within the very short 
time of the test. About half of the increase in strength is attributed 
to a dilatancy effect and half to the high rate of loading alone, 

It is a well known fact that when a dense sand is sheared it 
increases in volume. If shear occurs under conditions that prevent 
the movement of pore water, such as in the undrained triaxial test or 
in very rapid loading, the increase in volume will result in a decrease 
in pore pressure and an increase in strength, as shown by the equation 

s = (6 - u)tang 
in which the value of u (pore pressure) is decreased. These results 
are summarized in figure (19) and apply to fine and coarse sands for 
very short loading times. 

If the results of the tests made by Seed and Lundgren are 
applied in equation (15) the following results are obtained: 

Dynamic Resistance (Ry) = 1.40 x Static Resistance (Rg) 
or, Ry = Rg(1 + JV) 
where J is the damping factor and V the velocity of deformation, 


1.40 Re = Rg({l + JV) 


3 = 2.40-1.00 
= 


From figure (20), the load curve from one test, the velocity of 


deformation is estimated to be about 35 inches per second, and 
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substituting this value in the last equation gives 
is = a7 ios 
for the case of rapid loading in saturated sand. This value agrees 


closely with the value of 0.15 for J suggested by Smith (ref. 28). 


5. Work of Cunny and Sloan 


some preliminary tests were conducted by the Waterways Experi- 
ment Station, Corps of Engineers, Vicksburg, Mississippi in which small 
footings were dynamically loaded in a test tank filled with clay or 
dry sand (ref. 11). In these tests a loading machine capable of 
applying a load in a time of from 3 to 150 milliseconds was programmed 
to apply a preset load during a predetermined loading time. As a 
control, static tests were also performed on the same soil under 
identical conditions except for loading time. 

Dynamic and static tests were conducted on a sand compacted to 
a relative density of 96 per cent. Figure (21) shows the relationship 
between load and deformation for the static and one dynamic test on a 
J-inch-square footing. The load time curve for the dynamic test is 
shown in figure (22). From this curve it is seen that the loading 
time was 11 milliseconds and the maximum load reached was 3500 pounds. 
Since the maximum value of the load was limited by the setting of the 
machine, it is not necessarily the full strength of the soil under 
dynamic loading. 

By relating the values of static and dynamic loading at equal 
deformations, it is possible to arrive at an indication of the value 


of the damping factor under these conditions of loading. At a deformation 
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of 0.10 inch the static load is found to be 1600 pounds and the dynamic 
load 2650 pounds from the load-deformation curve of figure (21). The 
velocity of deformation at this point is estimated to be about 0 inches 
per second, or 3.33 feet per second. The corresponding value of J 
using equation (15) as before is found to be 0.197. 

6. Comparison of Dynamic Tests with Computer Solutions 

In order to determine whether the dynamic tests described above 
are qualitatively comparable with the equations used for computation 
of the total resistance, an analysis has been made of several sets of 
computer results in the same terms as the data obtained from the 
laboratory results, This analysis could only be carried to the point 
of maximum deformation, since the program was established to stop at 
this point. 

Figures (23) and (2h) are plots of the computer results for 
pile tip resistance-time, deformation-time, and load-deformation 
relationships for a 53-lb. H pile 0 feet long and a concrete pile 36 
inches in diameter and 176 feet long. The H pile is driven with a 
Vulcan #1 hammer and the concrete pile is driven with a Raymond )/0 
hammer. The tip resistance-time curve for the H pile, figure (2h), 
compares closely in shape with same curve for the dynamic laboratory 
test of the 9~inch-square footing shown in figure (22). This would 
tend to indicate that the load-time relationship resulting from 
equation (15) is a fair representation of the actual dynamic load 
mechanism, at least with respect to a dynamic test on a small footing. 

The pile tip velocity-time curves in figures (23) ad (2h) 


show a velocity of about 2 feet per second at the time of maximum 














— > & . 
ee J 
ee eee ee 
28 eT ce me one tet GL 
pl ™ " Gar eu os | oa 

4 1s 





eee Eh ee= —- 
ee it = eit & Ow 
= «ee as «© — ee! 
CF — = ee © en + a oe of Oh 
—MainW rE — se aa ee : 


4a =a 6 nm 4 «© 


“ia ele © -— £6 — za? 


Pte ee ee cen ime 
[S| i=. « - '«'@ @16' Gee ea 
> «- ~ whe * oe ee | 
= i ni Gale, Gat | 
SP pene! el ae a 


- > 
> © JP aa) - Se = ee oe 





SS —=_ 2 « *' «4 a, 
— —_— <= — * o le! / 
he meee“ — 
» « ° an on 


— — . - = & > = -—_—__> 





eae Saat ~~ ee ee 
aay Se ee tt 


aS y= Sfe-}=9 =| =e eels 
— Re ‘acing “tip -restet BD en ‘yelocity-anc plecanent—— ae ee = 
Seal: oe = iwith: tine for the: rah of 12 BP 53 pile| as computed by! the eee us 


—- 
ee = are oie ees, ee eee ae sfee 
—+1- re game Fe ee Se —— hb . t a fa ee o> + t { 
- : - —t tat 4 . me ee TH17 - =~? I> tie ah pe fetatal 
ae can oe eet ane pele wt Sieger —~y _— a eet a ——— woe b aoa 
' I 
—_ — o- -_ ce * ~ = = — ade © eee ia bmd—m! eo Jott = + + 4 Th 
= 5 = Son. » |= EB ne Tea . . z too pat tit ct cee. a apes ees erarat {ost 
tee | oh = eee ot eas > Beet eee 3 oo Hod 
soa ee = oe eet eee ~- af my ; cose o a fea fi per ++ a Sed 
a: eek oe so =. - = wd 5 ‘e « —! = * - ; ~ ~~ A — + vo eer iw -~++ — a m4 -+ he pa = oiadiandene aes 
Rn a Oe ho : tat = : ‘ a : $ : 1 : ete a 4 aeetoe ’ : 
: : ; : 1744 oe aos Talc es = , peor 
- a i” - - - - -- a : -_— 7 ~ |--— oo or -- = =4e-- + fp -wntg se f-- aaron = — 2 
. T . ' : 5 Ai ' ' 1 ’ . 4 T ‘ es , * z 1 ’ ° ? 
~ _ oo - 4-—-—-7---+4-- - ~ a , _~ —;>-c4e-e —- - - —} ++ “ . a 4b bay eet — —p - —-4~ - --+ ty et é — 
-* - a - §w- -~- + — : as ot fe - - i—— -§ = + fae t - we te + =: ae —++ 
1 . ' + (4 - ’ ‘ + * * sree * . } 
= eee fee a PS se Heel ae errata Senne op ~ +44 -* tis ta L$ - +4 nee ; as r ote any 
= : ; a cee eel Geb enk eels eeeceoneees er ntewien pant d. bo. ‘ . ~+— = : [ Be te ee eee 4 f+ f .—4— ; . 
te : I ¥ \ ; i aed i : ‘ TT. 1 
ri . ee! : 2 : ° =a t Tate wg 7, eee ala = an . ' Wary ' oo i? t } +++] T op i 
: T 1 ‘ia 1 ‘ 1: [AR 1 wary ? sales t it cy % ! 
: ae = shes any roieiaeat Llp ¢ . ; iinet re en eat 7 oe Ey } Sohbet : ie 
ar . Tee . ~ : Saar Tae ate - 71 — 3 a = —_ ———_— ~~ - -\- 
4 - r —_ . —_ = : ie cig mes 2 
ieee etal Saab aiaear gas is epee ant Lg ; aie Tt : : : a ad : 1 at 
: - = ’ . + +f? ~— ++ - = 7 
; ‘ ; P: ' ) : : 
cra eee aire rare a ’ : Print & 
’ y * Sa vatietie lear 
, ee + 
: - 
< _ 


pee yas aber a fee 
ees es ee ete a a eee padas Seb eee 
rm. 
Seapets HH anaes 
eee : 


HEE 

















i 
. = 
i— 4 = , nfl 
eee tteenw ds aes 
iets ae 
toi — + . 
erie aie ee 
1 ! . ei 
+ — 9 , fe 
= | ' i 
é + $e 
i 


E 


SaaS 
estas heey 
EEE 
falas 
SoH 


HEN 
re 


ENN 


\ 


~$4- 
fe 


1 —3- 9=10.15 § - 
==j-= es os 


a 

eo = —_ te — - —_— -_—— 
ae -O —-—-—- 9 —J— - + eer - 
BRT = == 


3 “Tt ea 7 


ENGR 
Ee 


2 
LEN, 
ay 


Ta 
THEE 
a 
mn 


a . 
foes 
ES ee 
{—-=- 5 





Sot tae 
Petes 
Eire Ee pee eee 
eet eee sess Eat 
sete PES Soe ee Co He 
aia Peeper SEE 
Se eas faneesiaed== eat 
ro aster ich Bete that Beate 
Es SS = pects pa= ee 
iS ve or 24 ia PERE eet Co 








60 


resistance for the 0-foot H pile and a velocity of about one foot 
per second for the 176-foot concrete pile. If it is accepted that 
most of the increase of dynamic resistance over the static resistance 
ls due to the increase in compressive strength of the soil due to 
rapid loading, and if this increase amounts to a factor of 1.0 for 
clays and 0. for sands, then the damping coefficient must have a 
value such that when multiplied by the velocity of deformation at 
the time of maximum load results in values of 1.0 for clays and 
0.4 for sand. Figures (23) and (2h) indicate that tip velocities at 
the time of maximum resistance are between 1,0 am 2.0 feet per second, 
and based on these values, the damping factor (J) would have a value 
ene0.5 to 1.0 for clay and 0.1 to 0.) for sand. 

If load vs. deformation is plotted from the information 
available in the curves of load and deformation vs. time for the 
12 BP 53 pile shown in figure (23), a curve shown in figure (25) is 
obtained. This curve is very similar in shape to the curve of dynamic 
load vs. deformation from rapid load tests performed on sarid in the 
laboratory as shown in figure (21). The fact that these curves are 
Similar in shape tends to indicate that the variation of load with 
deformation in the wave equation solution is of the same farm as 
that occurring in laboratory load tests. From this it is concluded, 
tentatively, that the computation of total driving resistance at the 
tip by the equation 

R= (D, - D,) x! (1 + Jv) 

is a valid representation of the actual physical process of rapid 


penetration of the soil by the pile. 
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7. Ground quake, which has been defined as the limit of elastic 
deformation of the soil, is normally thought of as a static concept, 
but it is an essential element in wave equation computations. The 
elastic ground compression occurring during static loading is not 
necessarjly the same that occurs during pile driving, and no attempt 
has been made to relate the two. Therefore, quake as used in this 
investigation is a dynamic property. 

A method of measuring quake during pile driving is described 
by Chellis (ref. 6). Essentially, this method measures the total 
elastic compression of the pile and soil. Pile compression alone is 


computed from the formula C, =R, L/AE, where R, is the ultimate 


2 

bearing load on the pile, Lis the length of the pile, A is the cross- 

sectional area, and © is the modulus of elasticity. This value is 

subtracted from the measured compression of pile and soil to obtain 

the value of quake. This method of measuring quake assumes that the 

pile is in compression simultaneously along its entire length, which 

is somewhat erroneous in view of the longitudinal stress wave concept. 
In extensive studies made of pile driving, Hiley made measure-~ 

ments of quake and reported that quake varies considerably for different 

conditions of ground and may average 0.1 inch in firm gravel, up to 

0.2 inch in firm clay and may have values as high as 0.5 inch in 

peaty material (ref. 17). With respect to driving conditions, Hiley 

summarized his observations as follows: 


Driving Conditions Quake Stress in Pile 
asy driving 0 - 0.10 500 psi 


O 
Medium driving 0.10 - 0.20 1000 psi 
Hard driving 0.10 - 0.30 1500 psi 
Very hard driving 0.0 - 0,20 2000 psi 
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It was also stated that quake may be as much as double these values 
if soft ground is immediately below the pile tip. These values 


for quake are used as a limiting range for investigations of load 


tests. 
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Chapter VI 


CORRELATION OF WAVE EQUATION RESULTS WITH PILE LOAD TESTS 


A, Methods of Approach and Effects of the Various Parameters 


Although there are many load tests published, there are com- 
paratively few which have not left out some vital part of the in- 
formation required for a thorough evaluation of the total pile driving 
and loading process. Partly due to this limitation, two general 
methods of approaching the problem of correlating wave equation 
computations with load test data were tried. 

The first method was particularly applicable to those tests 
for which very complete information was published. In these cases 
it was thought that most of the variables could be fairly accurately 
evaluated so that the investigation could be limited to a trial and 
error approach in which the principal unknowns could be varied over a 
wide range of values and tried in various combinations. Generally 
this means that ground quake and damping constants and side friction 
were varied with all other factors kept constant for the case in 
question. 

The second method of approach was adopted to try and derive 
as much information as possible from published load test data which 
were not reported as completely as those investigated by the first 
method of approach. It also considered desirable to run a series of 
data to study more fully the effects of the various parameters on the 
resistance-set relationship for particular piles and hammers. It 


was thought that both of these objectives could be met by selecting 
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@ particular pile type for which a number of reasonably good test 
data were available. Furthermore, it was desired that the pile meet 
the following requirements: 

(a) Be constructed of a material of which definite information 
was available as to its properties, 

(o) Be constructed of a material which of itself would least 
affect the soil around it and vice versa. 

(c) Have uniform properties throughout its length. 

(d) The pile as driven would be required to sustain the static 
loads to be imposed upon it. 

Tne reason for these requirements was to eliminate as many of 
the unknowns from the problem as possible. The steel H pile and steel 
pipe pile each satisfied these requirements better than other types. 
Because there were more usable load tests readily available for steel 
H piles, particularly the 12-inch H bearing pile weighing 53 pounds 
per foot, this shape was selected. It had the added advantage that 
most of the records available indicated that one type of hammer had 
been used, the Vulcan Number One; however, tests were available with 
other size hammers so that a further check of the method could be made, 

The main effort has been directed toward establishing correla- 
tion for each load test case within a limiting range of values of each 
of the variables involved. Correlation of wave equation computations 
with the actual pile load test involves the evaluation of five different 


variable factors. These are: 
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(a) Ground quake; 
(b>) Point damping factor; 
(c) Side damping factor; 


(d) Distribution of resistance between point and sides of 
the pile; 


(e) Per cent of static resistance acting during driving. 

These factors have been discussed in above sections with the exception 
of the last named variable, the per cent of static resistance acting 
during driving, which will be discussed below. 

In correlating theoretical calculations with load tests and 
the corresponding driving records, all variables except one must be 
held constant while varying one factor at a time. From practical 
considerations it was not possible to treat all five factors alike as 
variables, Since little is known of the nature of the side damping 
factor and sime the effect of side damping in most cases is relatively 
small compared to point damping, the side damping factorwas eliminated 
as a separate influence by assuming it equal to one third of the point 
damping factor for most cases. Smith indicated his belief that side 
damping would be less than point damping for a given soil, and he 
recommended a value which was one third of the point damping coeffi- 
cient (ref. 28). 

Although up to five different distributions of resi stance were 
investigated for each correlation case, this factor was assumed known 
when it came to the actual correlation using the first method. For the 
other method, distribution was treated within a range of values. The 


various distributions of resistance used for investigation are: 
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(a) End bearing 

(0) 75 per cent end bearing and 25 per cent side friction 

(c) 50 per cent end bearing and 50 per cent side friction 

(d) 25 per cent end bearing and 75 per cent side friction 

(e) All side friction 
Since the pile is generally divided into sections 5 to 10 feet long 
for purposes of calculation, the side friction acting on the lower 
section of the pile is indistinguishable from the point bearing resist- 
ance, except for damping, and adds to it for all practical purposes. 

Three methods of distributing side friction along the embedded 
length of the pile have been used. The method used most, and the 
simplest, assigns a uniform value of resistance for a unit length of 
pile regardless of depth. This method is termed rectangular distribu- 
tion. In another method the side friction is distritited in a manner 
so that the unit resistance is directly proportional to the depth, and 
is called triangular distribution. Both of these methods are included 
in the computer program, and selection of the applicable method is . 
made by use of a control statement in the program. Except for very 
long piles, the difference in effect on the resistance-set relationship 
of the two methods is small, so the first method was generally used. 
The third method of distributing frictional resistance is by estimating 
the actual resistance developed in each of the soil strata that the 
pile penetrates, 

The per cent of static resistance acting during driving is a 
factor that is influenced by time effects. In a sensitive clay, pile 


driving remolds the soil resulting in a reduction in strength during 
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and immediately after driving of the pile. If the load test is con- 
ducted after the soil has had time to regain its strength, the value 

of the static resistemce should be reduced when making the pile-driving 
analysis. Another type of time effect may result when a pile is driven 
into a sand layer through overlying layers of soft or medium clay. 

When the pile is test loaded, resistance initially developed in the 
clay may be transferred by plastic flow to the tip and lost from the 
standpoint of static load-carrying capacity. During driving, however, 
resistance in the upper clay layers will be present end will appear 

to result in heavy damping. Another effect is the liquefaction of 
saturated silts due to a quick condition being created by the driving 
action and resulting in greatly reduced resistance to penetration of 
the pile. The consideration of the effects of such factors in the 
calculations is very approximate and without complete soil data cannot 
be accounted for. 

In correlating wave equation computations with load tests in 
which sufficient data was available on which to base a reasonable 
estimate of the distribution of load between point bearing and side 
friction, the approach was to determine the various values of quake 
and damping that produced a resistance-set curve that included the 
point representing the ultimate load and final penetration of the test 
pile. Hypothetically, there would be an infinite number of such 
curves, but a limiting range of values was used for these variables 
that was determined on the basis of Smith's work and the experimental 
work on the dynamic properties of soil discussed above. The range of 


values for quake and damping investigated were: 
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(a) Quake 
(1) samt 0.05 to 0.20 
(2) clay 0.05 to 0.30 
(b) Damping 
(1) sand 0.10 to 0,20 
Cy clay OskO to 1,400 

These ranges of values were not rigidly adhered to, and in 
some instances outside values were investigated. The load test in- 
formation, input data used for the computer program for each load 
test, and results of the correlation attempts are included in tables 
F to BB. The accompanying figures are the ultimate resistance versus 
set (in blows per inch) curves drawn from the computer results and 
used as the basis for correlation. 

Similar curves used as the basis of correlation by the second 
method are included in appendix E, The length of piles used for the 
calculations was chosen to generally fit the load test data available, 
but since the effect of length was not clearly established the choice 
of 0-,80- md 120-foot piles was somewhat arbitrary. Calculations 
were made for 12 BP 53 piles in these lengths using a Vulcan Number One 
hammer with a hardwood capblock, a pile cap weighing 700 pounds, for 
all permutations of the ranges of ground quake, point and side damping 
factors and side friction distributions of a rectangular shape. For 
each combination of the variables, set in blows per inch was computed 
for each value of ground resistance ranging from hO to and including 
280 tons in 20-ton increments. The results were plotted so that on 


any graph there are five curves, one for each percentage of side 
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resistance assumed. For any one curve of ultimate resistance vs. set, 
all other variables are constant. 

These results show the following generalized effects of 
changing the variables indicated below: 

(a) Ground quake. Increasing ground quake increases the set 
in blows per inch, i.e. harder driving, for a given ultimate ground 
resistance for all values of side resistance including the point 
bearing case, 

(o) Increasing both the point and side damping increases the 
set in blows per inch for a given ultimate ground resistance for all 
values of side friction, including zero. 

{c) Side distribution--increases in per cent of side friction 
cause dramatic decreases in set in blows per inch for a given ultimate 
ground resistance, 

(d) Increased pile length--here the results are mixed. In the 
zero friction or end bearing case, the ultimate resistance tends to 
become the same for l0-, 80-, and 120-foot cases after about 2) blows 
per inch. Before that point is reached, the 80- and 120-foot piles 
have essentially the same curve of ultimate resistance vs. set, and 
the O-foot pile for a given resistance within this range tends to have 
an appreciably decreased set in blows per inch. For about 50% friction, 
the results may be practically identical for all three lengths. For 
higher values of friction the longer piles may have decreased set for 
given resistance. 

(e) Increased pile weight may result in decreased set in blows 


per inch for a given ultimate resistance; however, in the lower ranges 
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of resistance the difference may be slight. For example, for a side 
resistance of 25 per cent and all factors the same far two h5-foot 
H piles driven with a Vulcan O hammer, the difference in set for a 
14 BP 89 anc a 14 BP 117 is less than 2 blows per inch. The values 
are 8 and 6.6 respectively. 

(f) Increased hammer weight, but with the same stroke, results 
in decreased set in blows per inch for a given resistance and all 
other factors constant. 

(g) A softer capblock results generally in increased blows 
per inch for a given resistance, but the effect is really quite variable. 
Comparisons were made for an ultimate ground resistance of 100 tons for 
the following piles driven with a Vulcan Number One hammer: 

1) Concrete Piles--18 inches square with lengths 
of 40, 80, and 120 feet. It was assumed that head packing 
would be used, but this same soft material was assumed for 
the capblock also and the results compared with capblocks 
of hardwood and micarta. There was no appreciable differ- 
ence in sets for any of these conditions for any of the con- 
crete piles. 

2) Wood Piles--for the 80- and 120-foot lengths, the 
difference was slight for the three capblocks. For the 0- 
foot pile the difference between a micarta and a hardwood 
capblock were slight, but the very soft capblock made about 
10 blows per inch difference in the results. 

3) Steel piles--the capblock had essentially no effect 
on the 80- and 120-foot length piles, but the very soft cap- 
block made about 2 blows per inch difference with the 0-foot 
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The effects produced by varying some of these parameters are 
not surprising, but in some instances they may be, especially if one's 
judgment is based on results produced by popular dynamic formulas. 

For exanple, the results produced by the Hiley formula may indicate 

very slight differences between the cases when side friction is 

present and when it is totally absent. On the other hand, it is 
interesting to note that when side friction is a large percentage of 

the ultimate resistance to driving, the set vs. ultimate resistance 

curves may tend toward the curve produced by the Engineering News 
formula. The foregoing results produced by the wave equation method 

of analysis would seem to bear out its general suitability under widely 
varying conditions of driving equipment, pile type, and ground conditions. 

As can be seen by examining the following case studies, some 
of the data upon which to base correlation is rather complete, and 
some is not. The cases studied include various size piles, assorted 
pile driving equipment, in soils varying from gravel to clay, piles 
constructed of concrete and steel, different ratios of embedded length 
to driven length, and for varying combinations of soil strata, It is 
apparent that for many of the cases correlation can be obtained in 
several ways by various combinations of values for ground quake, point 
and side damping, and side friction. Because these factors are not now 
known with a high degree of certainty, the problem becomes one of 
searching for a pattern within the likely ranges of the variables for 
the various types of soils encountered, 

In some of the following cases it was necessary to estimate 


the values of Q and J which produce correlation, since the range of 
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values used in the computations did not extend far enough to include 
the load test point within computed load-set curves. The values for 
Q and J so estimated are based on the curves nearest to the load test. 
and the trend of variation between curves. It is believed that the 
estimated values are reasonably accurate considering the nature and 
accuracy of the basic data. 

Before examining the case studies in detail, the manner of 

identifying the set versus resistance curves is noted: 

Curve A -- 0% friction, i.e., an end bearing pile with no 
resistance to driving offered by the soil along 
the pile sides. 

Curve B -- 25% friction, i.e., 75% of resistance to driving 
offered by the soil at the pile point and 25% 
along the pile sides, 

Curve C -- 50% friction, i.e., 50% of resistance to driving 
offered by the soil at the pile point and 50% 
along the pile sides, 

Curve D -- 75% friction, i.e., 25% of resistance to driving 
offered by the soil at the pile point and 75% 
along the pile sides. 

Curve E -- 100% friction. All resistance to driving 


distributed along the pile sides. 
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TABLE F 
Case Number 2001 
Reference: (18) 
Pile: Prestressed concrete, 36 in. diameter, 176 ft. long 


Hammer : Raymond 3/0, 12500 1b. ram, 2 in. stroke 
Embedded length: hO ft. 

Oi): Soft to medium clay with point in stiff clay. 
Final penetration: 20 blows per inch 


Load test: Tested to failure at 190 tons including the weight of 
pile, tested one week after driving 


Input data: W(1) 12,500 Ram 

w(2) 8,800 (includes follower and driving cap) 

W(3)-w(20) 4,955 (based on 9.775 ft. segments) 

V 10.15 (ram velocity at impact for 2-ft. 
drop ) 

K(1) 10,600,000 (capblock spring constant) 

K(2) 2,180,000 (cushion and pile) 

ine) = (Falter 

K(19) 20, 750, 000 

K!(m) Soil spring constant depends on value of 


ultimate resistance which is varied. 


Distribution of resistance: 50% end bearing and 50% friction, based 
on computation of point capacity by the 
formula: P = 9AC. Friction considered 
to have a rectangular distribution. 


Results: See figures (27) and (28) 


Correlation is obtained for the following values 


of Q and J: 
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TABLE G 
Case Number 2002 
Reference: (18) 
Pile: Prestressed concrete, 36 in. diameter, 176 ft. long. 


Embedded length: 0 ft. 
Hammer: Raymond 4/0, 15,000 lb. ram, 34 in. stroke. 
Soil: soft to medium clay with point in stiff clay. 
Final penetration: Refusal, last 3/4 in. at 99 blows. 
Load test: Tested to failure at 230 tons including the weight of pile. 
Input data: W(1) 15,000 
W(2) 8, 800 
V 12.08 
(Rest of input data same as for case 2001) 
Distribution of resistance: Estimated between 75% and 50% point bearing. 
Results: See figures (29) and (30). 


Correlation is obtained for the following values 
of Q and J; 


Q J J! 
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TABLE H 
Case Number 2003 
Reference: (18) 
Pile: Prestressed concrete, 36 in. diameter, 176 ft. long. 


Embedded length: 0 ft. 
Hammer : Raymond 3/0, 12,500 1b, ram, 39 in. stroke 
Soil: Soft to medium clay with point in stiff clay 
Final penetration: 10 blows per inch 
Load test: Tested to failure at 165 tons including weight of pile. 
Input data: W(1) 29500 
V 269 
(Rest of input data same as for case 2001) 
Distribution of resistance: Estimate 50% end bearing and 50% friction 
Results: See figures (31) through (3h). 


Correlation obtained for the following values 


of Q and J; 
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TABLE I 
Case Number 200) 
Reference: (18) 
Pile: Prestressed concrete, 36 in. diameter, 80 ft. long. 
Embedded length: 3.5 ft. 
Hammer : Raymond 3/0, 12,500 1b. ram, 39 in. stroke. 
Soil: Medium and stiff clay with point in very stiff clay. 


Final penetration: 20 blows per inch. 


Load test: Tested to failure at 365 tons. 


Input data: W(1) L285 OO 
w(2) 8, 800 
W(3) - W(10) 5,070 
V 267 
K(1) 10,600, 000 
K(2) 1, 020, 000 
K(3) - K(9) 20, 300, 000 
e7 On6 
eo 0.5 


Distribution of resistance: From computations based on load test and 
soil data, it is estimated that 87% of the load was 
carried by the pile point at failure, Correlation should 
be made for distributions between 100 per cent and 75 per 
cent of the load carried by the point. 


Results: See figures (35) through (37) 


Correlation is obtained for the following values 


of Q@ and J: 
2 ones ha 
Om. 0.45 0.15 
.05 O26 OF2 


Remarks: In this case, a hole 2) ft. in depth and 42 in. diameter 
and below that for a depth of 27 ft., a hole 35 in, 
diameter was pre-excavated for the pile. The effect of 
the 35 in. pre-excavation is uncertain, but appears to 
have considerably reduced the dynamic resistance to driving. 
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TABLE J 
Case Number 2005 
Reference: (18), Test 7-I 
Pile: Prestressed concrete, 36 in. diameter, 68 ft. long. 


Hammer: Raymond 3/0, 12,500 1b. ram, 39 in. stroke 


Embedded length: 31 ft. 

Soil: otiff clay. 

Final penetration: lO blows per inch. 

Load test: Tested to failure at 260 tons. 

Input data: Same as case 200), except: 
W(3) = W(9) 4,930 
K(3) - K(8) 20,900, 000 


Distribution of resistance: 


From calculations based on load test and soil data it is 
estimated that 43% of the load was carried by the point 
of the pile failure. For correlation, 50% of the load is 
considered to be carried by the point. Friction is con- 
Sidered to be distributed over the lower 31 ft. of the pile. 


Results: See figure (38). 


Correlation is obtained for the following values 
of Q and J: 


Q J J! 
OFL0 ae OPES. 
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TABLE K 
Case Number 6006 
Reference: (1) 
Pile: Prestressed concrete, 5) in. diameter, 96 ft. long. 
Hammer : Raymond 3/0, 12,500 lb. ram, 39 in. stroke, wood capblock. 
Embedded length: 60 ft. 
Soil: Soft, medium and stiff clay and silty clay with some sand 


with pile of point in dense, fine, gray sand. See figure 
(39) for soil profile. 


Final penetration: 5 blows last 3/l in. (60 blows per inch). 


Load Test: Tested to 20 tons, failure had not occurred, Bayliss cal- 
culated a minimum ultimate load capacity of 98 tons based 
on load test data and soil characteristics using static 
formulas (ref. 1). 


Input data: W(1) 12,500 
W(2) 6, 000 
w(3) = w(12) 6,200 
V 12.9 
K(1) 2,715,000 
K(2) 2,025, 000 
K(3) - K(11) 27, 200, 000 
e7 0.50 
&4 0,50 | 
i 0.00033 


Distribution of resistance: 
Based on Bayliss! computations (ref. 1) it is estimated 


that 50 per cent of the load was carried as em bearing 
with friction distributed over the lower 60 ft. of the 


Dake. 
Results: See figures (39) through (2). 


Correlation is obtained for the following values 


of Q and J: 
ay a nee 
Onl 0.15 0.05 
0.15 0.16 0.05 


0.05 2: .05 extrapolated 
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TABLE L 


Case Number 6007 


Reference : 


Pile: 


Hammer : 


(19) 


19 


Pipe pile, 18 in. diameter, 3/8 in. wall thickness, 80ft. 


long. 


Vulcan No. 1, 5,000 1b. ram, 36 in. stroke. 


Embedded length: 


Soil: 


Final penetration: 


Load test: 


Input data 


Distribution of resistance: 


Results: 


Predominantly plastic clays to a depth of 80 ft. underlain 
by uniformly graded fine sand having a relative density of 


approximately 80 per cent. 
2 blows per inch. 


Tested to failure at 81 tons. 


Pile tip stopped 5 ft. above sand. 


: W{1) 5000 
W(2) 1,000 
W(3) - w(12) 562 
V 12.4 
K(1) 2,715, 000 
K(2) - K(11) 6,520, 000 
e] O70 
eo eo 
T 0.00025 


See figure (3). 


Load carried as friction. 


Correlation is obtained for the following values 


of Q and J: 
Q J 


Oe 0.15 


a! 
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TABLE M 

Case Number 6008 
Reference: (19) 
Pile: Steel pipe pile, 18 in. diameter, 3/8 in wall thickness, 

96 ft. long. 
Hammer : Vulcan No. 1, 5,000 lb. ram, 36 in. stroke. 
Embedded length: 75 ft. 
Soil: Plastic clays to a depth of 80 ft. underlain by a uniformly 


graded fine sand having a relative density of approximately 
80 per cent. The pile tip was driven about 5 ft. into the 


sand layer. 


Final penetration: 226 blows for the last 6 inches (37.7 blows/in. ) 


Load test: Tested to failure at 2lh tons. 


Input data: W(1) 5,000 
W(2) 1, 000 
W(3) - W(12) 676 
V 12.4 
K(1) 2,715,000 
K(2) - K(11) 5,420,000 
al 0.50 
eo Le 
T 0.00025 


Distribution of resistance: 


Based on two different load tests, one with the pile tip 
in sand and the other (case 6007) with the tip 5 feet 
above the sand layer, it was determined that 159 tons 
or 65 per cent of the load was carried by the point (the 


lower 5 ft.). It is assumed the frictional load is 


distributed over the embedded length, 


Results: See figure (hh). 


Correlation is obtained with the following values 


of Q and J: 


Q J J! 


0.05 0,10 0.03 
O70 0,05 0.05 (estimated) 
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TABLE N 
Case Number 6011 
Reference: (20) 
Pile: ly-inch H pile weighing 73 1b./ft. and 81 ft. long. 
Hammer : Vulcan OR, 9,300 lb, ram, 39 in. stroke, 
Embedded length: 80 feet. 
oi]: Alternating strata of silts, sandy silts, and silty sands 


with some interspersed clay strata for a total thickmess of 
about 50 to 60 feet. Clean sands from 0 to 60 feet thick 
lie beneath the silts, The pile penetrated 32 feet into 
the sand layer. 


Final penetration: 28 blows per foot (2.3 blows/inch) 


Load test: The pile buckled under a test load of 3h tons; at this 
point the pile had reached a gross settlement of 1.2 inches. 
From the load-settlement diagram it appears that the pile 
Capacity had been reached, based on the usual criteria, 


Input data: w(1) 9, 300 
, W(2) 700 
W(3) - W120 64,8 
V 12.9 
K(1) 2,715,000 
K(2) ~ K(11) 7,270, 000 
ey 0.50 
€9 ian 
iL 0,00025 


Distribution of resistance; 


This pile was equipped with strain rods in order that the 
distribution of load on the pile could be determined during 
the load test. The high value of set (low value of blows 
per inch) for this pile tends to indicate the possibility 
of liquefaction in the silt layers. Assuming this to be 
the case, correlation is made on the basis of the load 
carried to the sand layer which was 170 tons. 


Results: See Figures (5) and (6). 


Correlation is obtained for the following values 
of Q and J: 
Q J J! 


eal 15 .05 
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Case: 
pource s 
Pile: 
Embedment : 
Hammer : 


Soil: 


oi 
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TABLE 0 


Page 31, reference (16). 


65' 12BP53 

55. feet. 

Vulcan #1 

10' silty sand 

1k! loam 

26! fine sand, silty sand, and clay 
10' fine to medium sand 

10! fine to medium sand and gravel 


then fine sand 


Final penetration: 1,63 blows per inch. 


Load test results: 5 tons 17 days after driving. 


Results of calculations by the wave equation method: 


a. 


Comments: 


Lf Q 


G,10 
CAL 


it It 


J'= 0,03 


correlation is obtained with a side friction of 2%, 


If Q 


0.10 
0.20 


Mt 


Jt= 0,067 


correlation is obtained with a side friction of 17%. 


If Q 


0.20 
OG 
0.03 


correlation is obtained with a side friction of 2%. 


If Q 


O720 
0,20 


i it 


J'= 0,067 


correlation is obtained with a sice friction of 32%. 


Hiley Formula gives ultimate resistance of 52 tons. 
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TABLE DD 
3 
Page 34, reference (16). 
12BP53 
47 feet. 
Vulcan #1 


35' silty clay 
then sand 


Final penetration: 4.57 blows per inch. 


Load test results: 100 tons, but ultimate resistance not reached, 


Results of calculations by the wave equation method: 


Ae 


Comments: 


it Q = O,10 
@ = 0.10 
J' = 0,03 


correlation cannot be obtained for these values for a 
failure load of 100 tons; however, these are the results: 


22 tons for no, side friction 
146 tons for 25% side friction 
185 tons for 50% side friction 


If Q = O,10 

JI = 0.20 

J! = 0,067 
correlation is obtained with no side friction, 
If Q = 0.20 

J = 0.10 

J' = 0,03 
correlation is obtained with no side friction. 
If Q = 0,20 

J = 0.20 

J = (07067 


correlation is obtained with a side friction of 25%, 


It is not known just what the ultimate resistance would have 
been if the load test had been carried all the way to failure. 


The driven length of the pile was not given; 57' was assumed. 


Hiley Formula gives ultimate resistance of 96 tons. 
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Case: 
Source $ 
Pile: 
Fmbedme nt ; 
Hammer : 


Soil: 


3h 


TABLE P 
35 
Page 35, reference (16). 
12BP53 
Su.) feet. 
Vulcan #1 
15' sandy, silty, clay and loam 
10! sand 
10! sandy loam 
2)! sand 


then sandy loam 


Final penetration: 4.57 blows per inch. 


Load test results: 100 tons, but ultimate resistance not reached, 


Results of calculations by the wave equation method: 


Qe 


Hag; = 0,10 
J = Q.10 
J'§ = 0,03 


correlation cannot be obtained for these values for a failure 
load of 100 tons; however, these are the results: 


113 tons for no side friction 
hi tons for 25% side friction 
182 tons for 50% side friction 


ties = O.2.0 
J = Q,20 
J' = 0,067 


correlation is obtained with a side friction of 8% for 
a 100-ton ultimate resistance. Other results are: 


11h tons for 25% side friction 
151 tons for 50% side friction 


if Q = 0.20 
J = 0.10 
gt 0,03 


correlation is obtained with a side friction of 2% for 
a 100-ton ultimate resistance. Other results are: 


122 tons for a 25% side friction 
163 tons for a 50% side friction 
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Case 35 (continued): 


ree = 0,20 
c= O220 
J§ = 0,067 


correlation is obtained with a side friction of 26% for a 
100—ton ultimate resistance. Other results are: 


103 tons at 25% side friction 
129 tons at 50% side friction 


Comments: It is not known just what the ultimate resistance would 
have been if the load test had been carried all the way 
to failure. 


The driven length of the pile was not given; 65!' was 
assumed, 


Hiley Formula gives ultimate resistance of 92 tons. 
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Case: 
Source: 
Fide: 
Embedment : 
Hammer : 


Soil: 
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TABLE Q 


36 

Page 36, reference (16) 
12BP53 

49.8 feet. 

Vulcan #1 


15' sandy, silty, clay and loam 
then sand 


Final penetration: .15 blows per inch. 


Load test results: 100 tons, but ultimate resistance not reached. 


Results of calculations by the wave equation method: 


ae 


If Q = Q,10 
J = OQO,10 
g' = 0,03 


correlation cannot be obtained for these values for a 
failure jJoad of 100 tons. Other results are: 


108 tons for no side friction 
135 tons for 25% side friction 
178 tons for 50% side friction 


If Q = 0,10 
J = 0.20 
J' = 0,067 


correlation is obtained with a side friction of 7%for 
a 100-ton ultimate resistance, Other results are: 


115 tons for 25% side friction 
143 tons for 50% side friction 


If Q = 0,20 
J = Q.10 
J'§ = 0,03 


correlation is obtained with a side friction of 7% for 
a 100-ton ultimate resistance. Other values are: 


115 tons for 25% side friction 
152 tons for 50% side friction 
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Case 36 (continued): 


d. 


Comments: 


moe) - 0.20 
J = Q.20 
Jt = 0,067 


correlation is obtained with a side friction of 32% fer 
a 100-ton ultimate resistance. Other values are: 


91 tons at 25% side friction 
121 tons at 50% side friction 


It is not known just what the ultimate resistance would 
have been if the load test had been carried all the way to 
failure, 


The driven length of the pile was not given; 60' was 
assumed, 


Hiley Formula gives ultimate resistance of 92 tons. 
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Case: 
source > 
Pile: 
Embedme nt : 
Hammer : 


eel: 
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TABLE R 
67 
Page 67, reference (16) 
114' 1OHh2 
95 feet 
Steam, double acting, 950 pounds with stroke of 16 inches. 


Silt and fine sand, pile hit hard strata 


Final penetration: 32 blows per inch 


Load test results: 60 tons 


Results of Calculations by the wave equation: 


a. 


Comments: 


If Q = 0,10 
J = 0.15 
Jt = 0,05 


correlation cannot be obtained; however, 100 tons was the 
calculated value for no side friction, 


Siew. = 0.10 
J = 0,20 
J' = 0,067 


correlation cannot be obtained; however, 88 tons was the 
calculated value for no side friction. 


If Q = 0.20 

Jt = 0,05 

d = 0.15 
correlation is obtained with a side friction of 15%. 
If Q = 0,20 

J = QO,20 

Hom 007 


correlation is obtained with a side friction of 04. 
Time between driving and load testing is unknown. 


Hiley Formula gives an ultimate resistance of 25 tons, 
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Case: 
source : 
Pile: 
Embedme nt : 
Hammer : 


Boll: 
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TABLE S 


fa 

Page 71, reference (16) 

20! 10Hh2 

19 feet 

Gravity, 3,150-pound ram and 5-foot drop 


Sand and Gravel 


Final penetration: 2 blows per inch 


Load test results: 6h tons 


Results of Calculations by the Wave Equation: 


a. 


Comments: 


oo = 0,10 
Jee] 6.15 
Jt = 0.05 


correlation is obtained at a side friction of both 75% 
and 100%. 


If Q = 0,10 

J = 0,20 

J' = 0,067 
correlation is obtained for a side friction of 100%. 
If Q = 0,20 

t=) 0.15 

fe O05 
correlation is obtained for a side friction of 100%. 
If Q = 0.20 

J = 0,20 

J' = 0,067 


correlation cannot be obtained; however, the ultimate 
resistance is calculated to be 60 tons for 100% friction, 
or O tons for the end bearing case. 

Time between driving and load testing is unknown. 


Hiley equation gives an ultimate resistance of 80 tons. 
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Case: 
Source : 
Pade : 
Embedment : 
Hammer : 


Soil: 


90 


TABLE T 
e 
Page 75, reference (16) 
30! 1OH2 
23 feet 


Drop--3000-pound ram, 5-foot drop 


Fine Sand 


Final Penetration: 5 blows per 2 inches 


Load Test Results: 90 tons 


Results of Calculations by the Wave Equation: 


Be 


hee = O.10 
J = 0,15 
J§ = 0,05 
correlation is obtained for 100% side friction. 
If Q; = 0.10 
ioe] “0,20 
was 90,067 
correlation cannot be obtained for a failure load of 
90 tons; however, these are the results: 
5S tons for no side friction 
6h, tons for 50% side friction 
70 tons for 100% side friction 
few = 0,20 
io 
J' = 0,05 
correlation cannot be obtained for a failure Load of 
90 tons; however, these are the results: 
Sh tons for no side friction 
60 tons for 50% side friction 
60 tons for 100% side friction 
If Q = 0,20 
we= O720 
eee OOO] 


correlation cannot be obtained for a failure load of 
90 tons; however, these are the results: 


7. rer eae 
eS 
-_ im 
— ee | ee 


mes te een 2 Uae 7 


| 
=i mt <t © emtnetinces® € 
** gy "a 


> a0 





i °* _ "— "ot te 
‘7 i . 2) See) 
oun ke wey earn 
“", ae ad 
a r 
: . = = a e tl pe 


m0 —_——e iammed to 








Comments: 


TABLE T (cont'd, ) 


48 tons for no side friction 
55 tons for 50% side friction 
60 tons for 100% side friction 
Time between driving and load testing unknown, 


Hiley formula gives Ru of 80 tons. 
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TABLE U 
Case: 76 
Source: Page 76, reference (16) 
Pile: hO' 1OH2 
Embedment: 3h feet 
Hammer : Drop--3000-pound ram, 5-foot drop 


Soil: Fine Sand 
Final Penetration: 2 blows per inch 
Load Test Results: 75 tons on second loading, 50 tons first loading 


Results of Calculations by the Wave Equation: 


a. If Q = O.10 
wee= 0.15 
J' = 0.05 
correlation is obtained with side friction of 25%. 
Eee ir © = Q;10 
J = O,20 
Jt‘ = 0,067 


correlation is obtained for a side friction of 35% (the 
difference in this case between 25% and 50% friction 
is only 1 ton). 


Cott O = O,20 
J = 0,15 
Jt = 0,05 
correlation is obtained with a side friction of 75%. 
ad. If Q = 0,20 
J = 0,20 
ies 0,067 


correlation is obtained with a side friction of 90%. 


Comments: Time between driving and load testing unknown, 


Hiley formula gives an ultimate resistace of 60 tons. 
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TABLE V 
Case: 89 
Source: Reference (12) 
Raikes 4S' J)BP89 
Embedment: lj)! 
Hammer: Vulcan 0 
Soil: Same as case # (1) 
Final Penetration: 30 blows per inch 


Results of Load Test: 300 tons 


Results of Calculations by the Wave Equation Method: 


Comment: 


Correlation is obtained with 25% friction and 


Q = 0.20 
J = 0.20 
Jt = 0,067 


This pile was tested 5 days after driving. The 300~ton 
load appears to be close to the ultimate resistance of 
the ground, Load increments were 75 tons, 150 tons, 
225 tons, 225 tons (repeated), and 300 tons. 


Hiley formula gives ultimate resistance of 230 tons. 
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TABLE W 
Cases tet 
Source: Page 113, reference (16) 
Pile: 4O' precast c oncrete weighing 2.7 tons 


Embedment: 31.5! 
Hammer : Vulcan #1 
Driving Cap: 0,525 tons 


Soil: 18! coarse sand 
26' fine, yellow sand 


Final Penetration: 10 blows per 2.05 inches 
Load Test Results: 100 tons 


Results of Calculations by the Wave Equation: 


eet 2 = O,10 
= eS 
J! = 0,05 
correlation is obtained for 50% friction. 
Celio. = O,10 
J = 0,20 
J! = 0,067 
correlation is obtained with a side friction of 85%, 
Coat GO = 0,20 
J = 0.15 
t= 0.05 
correlation is obtained with a side friction of 75%. 
d, if 2 = 0.20 
J = O.20 
J' = 0,067 


correlation is obtained with a side friction of 95%. 
Comments: Time between driving and load testing is unknowm, 


Results using the Hiley equation are 122 tons. 
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Case: 
source: 
Pile: 
Embedment : 
Hammer ; 


Soil: 


95 


TABLE X 


117 
Reference (12) 
4St 14BP117 
lh feet 
Vulcan O 
9' sand & gravel 
12! fine to medium sand 
15! sand and gravel 
6! sandy silt 


2' weathered, silty shale 
then hard, silty shale 


Final Penetration: 22 blows per inch 


Results of Load Test: 300 tons 


Results of Calculations by the Wave Equation Method: 


Comments: 


Correlation is obtained with 25% friction and 


Q = 0,20 
dae 0520 
J' = 0,067 


This pile was test loaded 2 days after driving. The 
300-ton load appears to be close to the ultimate 
resistance of the ground. Load increments were 75 tons, 
150 tons, 225 tons, 225 tons (repeated), md 300 tons, 


Hiley formula gives ultimate resistance of 10 tons. 
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Case: 
oource ; 
Pile: 
Embedment : 
Hammer : 


Soil: 


96 


TABLE Y 


Lhd 

Page 14h, reference (16) 
120! 12BP53 

61.5! 

Vulcan #1 


5! clay loam 
ll' sandy, silty, clay 
12! medium sand 

5' coarse sand 

39' medium clay 

10' stiff clay 

10! very stiff clay 

10' medium clay 

ll! hard clay with pebbles 
then medium clay 


Final Penetration: 8.0 blows per inch 


Load Test Results: 117 tons ultimate resistance 


Results of Calculations by the Wave Equation Method: 


Be 


Comments: 


If Q = 0.10 
ee Oe 0 
J'§ = 0,03 


correlation cannot be obtained for a failure load of 117 
tons, At full point bearing (0% side friction) the 
ultimate resistance is calculated to be 132 tons, 


ioe ©: = O220 

J = 0,20 

J' = 0,067 
correlation is obtained with a side friction of 5%. 
fH. 8 = 0.20 

J = O,10 

Jt = 0,03 
correlation is obtained with a side friction of Of. 
If Q = 0,20 

JI = 0,20 

OOO | 


correlation is obtained with a side friction of 20%. 


Hiley Formula gives ultimate resistance of 90 tons. 
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TABLE Z 
Case: 15 
Source: Page 145, reference (16) 
Pile: 106! 12BP53 
Embedment: 105! 
Hammer Vulcan #1 
Soil: Same as case # lh) 


Final Penetration: 10,3 blows per inch 
Load Test Results: 182 tons ultimate resistance 


Results of Calculations by the Wave Equation Method: 


a. If Q = 0,10 
J = O,10 
J' = 0,03 


correlation is obtaired with a side friction of 19%. 


b. If Q@ = 0.10 


J = QO,20 

J' = 0,067 
correlation is obtained with a side friction of 35%. 

c. If @@..= ‘Or20 

J = 0,10 

J' = 0,03 
correlation is obtained with a side friction of 36%. 

ad. LiaeQ =4%0.20 

J = 0.20 

J' = 0,067 


correlation is obtained with a side friction of 7%. 


Comment: Hiley Formula gives ultimate resistance of 86 tons. 
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TABLE AA 
Case: 1h6 
Source: Page 146, reference (16) 
Pile: 100! 12BP53 
Embedment: 97.1! 
Hammer : Vulcan #1 
Bod]: Same as Case #1) 


Final Penetration: h.75 blows per inch 
Load Test Results: 91 tons 
Results of Calculations by the Wave Equation Method: 


ae If Q = 0.10 
J = 0,10 
J§' = 0,03 
correlation cannot be obtained for a failure load of 
91 tons. At full point bearing (0% side friction) the 
ultimate resistance is calculated to be 110 tons. 


poli. = OOO 
J = 0,20 
J' = 0,067 
correlation is obtained with a side friction of 2%. 
ec. If Q = 0,20 
J = 0,10 
J'§ = 0,03 
correlation is obtained at 2% side friction. 
d. If Q = 0.20 
J = 0.20 
J! =90 067 


correlation is obtained with a side friction of 17%. 


Comment: Hiley Formula gives ultimate resistance of 72 tons. 
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Case: 
Source : 
Pile: 
Embedment : 
Hammer : 


Soil: 


a9 


TABLE BB 


ier 

Page 161, reference (16). 
50! 12BP53 

32 feet 

Vulcan #1 


geet i 
25' coarse sand 
18' medium sand 
15! fine smd 
9' medium clay 
20' very stiff clay 
22' medium clay 
then very stiff clay 


Final Penetration: 4.83 blows per inch. 


Load Test Results: 167 tons 


Results of Calculations by the Wave Equation Method: 


ae 


Comment: 


If Q = 0.10 

J = 0,10 

J‘ = 0,03 
correlation is obtained with a side friction of 35%, 
If Q = 0.10 

J = 0,20 

Jt = 0,067 
correlation is obtained with a side friction of 55%. 
If Q = 0.20 

J = 0,10 

J' = 0,03 
correlation is obtained with a side friction of 53%. 
If Q = 0,20 

J = 0.20 

Jt = 0,067 


correlation is obtained with a side friction of 70%. 


Hiley Formula gives ultimate resistance of 96 tons. 
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C, Discussion of Results. 

Each of the cases correlated by the first method are dis- 
cussed below: 

1. Case number 2001, table (F). This and the next two cases 
are 36-inch diameter prestressed concrete piles driven at Lake Maricaibo, 
Venezuela. The load carrying capacity of the pile tip was calculated 
using the equation, P = 9AC, and a value of cohesion (C) of 2560 p.s.f., 
giving a capacity of 163,000 pounds. This value is 2 per cent of the 
total load carried by the pile at failure, and the assumption that the 
load is carried one half in friction and one half by point bearing 
would be reasonable, If the load carried as friction is considered to 
be evenly distributed on the embedded length of the pile, each ten-foot 
section would develop a resistance equal to 47,500 pounds. If the 
distribution of the side friction resistance were assumed to vary 
directly as the depth, the distribution would be as shown in figure 
(26), which also compared both distributions with that estimated from 
the shear strength of the soil. In the calculated values the full 
value of the shear strength of the soft clay layers was assumed to be 
effective and a reduced value of 800 p.s.f. was used for the stiff 
clay layer between elevations -153 and -143 ft. This value is con- 
sistent with Nordlund's estimate (ref, 23) and the findings of others 
that have investigated the friction valves of piles in stiff clay 
(ref. 33, 38). Comparing the load distribution in figure (26), it 
appears that the rectangular distribution is as good a representation 
of the estimated actual load distribution as the triangular method 


would be. From figure (27) it is seen that values of Q of 0.1 and 
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J of 0.9, and Q of 0,2 and J of 0.3 produce nearly identical curves 
that give a resistance value slightly higher (200 tons vs. 190 tons) 
for the set of 20 blows per inch. 

Case 2002, This is the same pile as case 2001 redriven to 
refusal with a heavier hammer, The load capacity was not increased 
a great deal by the redriving, going only fram 190 tons to 230 tons. 
Correlation was obtained only with a quake of 0.20 using a value of 
J of 0.80 as shown in figure (29). This is quite a departure from the 
value of J obtained in the previous case and would tend to indicate 
that for harder driving, quake increased. Correlation could not be 
obtained for a quake of O.1 with values of damping as high as 1.2. 

Case 200i. The tip of the pile in this case penetrated into 
clay having a shear strength of 10,000 p.s.f. as determined from 
unconfined compression tests. Results of carrelation ted to indicate 
a smaller value of quake than experiencedin the previous cases. 

Figure (38) shows the results from wave equation computations. 

Case 2003. In this case correlation was obtained with values 
of 2 and J of the same approximate values as in case 2001 as indicated 
in figures (34) and (35). Figures (31) to (33) are included as a matter 
of interest, and show that correlation could be obtained with different 
distributions of resistance, 

Case 2005, This pile is of the same type as the previous one 
but was driven into a somewhat softer clay having a shear strength of 
about 3000 p.s.f. Figure (39) is the basis of correlation, and indicates 


values of J in the range more-or-less expected. 
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Case 6006. The pile in this case was not loaded to failure, 
but, from analysis of the load test and driving records, the engineer 
Supervising this test estimated the value of side friction and point 
bearing resistance developed under static loading. From this the 
ultimate load carrying capacity of the pile was computed to be a mini- 
mum of 498 tons. This value was used in correlating the load test 
with wave equation computations. The correlation for a value of Q of 
0.20 was obtained by extrapolation, but it is believed sufficiently 
accurate considering the general order of accuracy of the other vari- 
ables and data. Figure (0) is the load settlement curve obtained 
from the load test, and figures (1) and (2) show the basis for 
correlation. 

Cases 6007 and 6008, These are two pipe piles driven through 
80 feet of plastic clay underlain by a thick sand stratum. One pile 
was driven about 5 feet into the sand and the other was stopped 5 feet 
above the sand layer. The object of the load tests was to determine 
the load carried by the point in sand. Correlation for the pile that 
extends into the sand (case 6008) was obtained in the lower end of the 
expected range of values for the point damping factor and quake, with 
values of 0.05 and 0.10 respectively. No attempt was made to account 
for the fact that the pile itself is largely embedded in clay. For 
the pile stopped above the sand layer, correlation was obtained for 
values of Q of 0,10, J of 0,15 and J! of 0,05. Since the load is 
carried entirely by friction, all of the damping is due to J‘, the 
side damping factor. These values are not consistent with the expected 


range of values for clay, and indicate that the value of quake is smaller 
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than expected, whigh would be consistent with the easy driving condi- 
tions, that the relationship between J and J' is in error, or that 
side damping does not develop when a steel pipe pile is driven through 
plastic clay. With respect to case 6003, pile point in sand, it is 
possible that the clay was somewhat sensitive and regained strength 
during the period before the load test was made. If the sensitivity 
of the clay were such that one half of its strength were lost dw to 
remolding and then regained after a period of time, the ultimate re- 
sistance during driving would be reduced by about O tons. This would 
produce correlation at avalue of Q of 0.10 and J of 0.15. 

Case 6011. The high load capacity developed by this pile and 
the low resistance to penetration indicate that a quick condition:in 
the overlying silt layers occurred. Ifthis were the case, the resistance 
to driving would be developed on the portion of the pile embedded in 
sand with no resistance occurring in the silt layers. Correlation is 
obtained on this basis as shown in figures (45) and (6). 

The results of correlation attempts for these cases are sum- 
marized in figure (67). The piles driven with the tip in sand are 
closely grouped and fall within a range of values for J of 0.05 ta 0.2, 
and for values of Q of 0.05 to 0.10. This result tends to confirm 
Smith!s recommended value for J of 0.15, for piles driven in sand. 

Results of correlation for piles driven in clay are spread 
over a range of values for J of 0.15 to about 1.4, and for values of 
Q of 0.1 and 0.2. Case 6007, which has a J value for correlation of 
0.15, is a 100% friction pile, and the point damping factor (J) has no 


influence. This possibly indicates that the one-third relationship 
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Blows} % Load Values of J for 
per at point | Q= 0.05 Q= 0.1 


J 





Test 
Load 





190 | 20 50 0.9 
132 50 = 75 
rome) 50 1.0 
ey S025 0.45 
40 | 100 = 75 5.0 
60 | 50 0.2 0.15 
ss 0 0.15 
38 75 = 50 op) 0.05 
eS | 200 ="75 0.15 





0.4 0.6 0.8 1.0 
Point Damping (J) 





Summary of Correlation Results - First Method 


Figure (67) 


10h, 


between J and J! is too gross an assumption in cases where all or most 
of the load is carried by friction. Im general, the results in clay 
are consistent with values of J expected, based on dynamic tests made 


on soils in the laboratory, 


D, Summary of Correlation Results for Second Method of Appr oach 


For the cases studied, the following values of ground quake, 
point damping, and side damping seem to be appropriate for the soil 
types indicated. The percentages of side friction given are for steel 
H piles unless otherwise indicated. Furthermore, the principal case 
studies used in arriving at the values are shown: 


1. Coarse Sand: (Case Numbers 113 and 161) 


Q = 0,10 
J = 0,15 
Jt = 0,05 


Side Friction is 35%. 


2, Sand and Gravel Mixed: (Case Number 71) 


Q = 0.10 
J = 0,15 
Jt = 0,05 


Side Firction is between 15% and 100%. 


3. Fine Sand: (Case numbers 75, 76, and 113) 


Q = 0.15 
2. = 0.15 
yt = 0,05 


Side Friction is 100% (use 50% for straight-sided 
concrete pile). 


4. Sand layers combined with layers or strata of clay or 
loam or both, but sand layers predominate: (Case numbers 


31, 34, 35, 36, lk, 145, 146, and 161 apply) 


Q = 0,20 
J = 0.20 
J! = 0.067 


Side Friction is 25% 
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5. Silt and Fine Sand underlain by Hard Strata: (Case 


number 67) 
Q = 0,20 
J = 0,20 
Jt = 0.067 


Side Friction is 04. 


6. Sand and Gravel Underlain by Firm Strata such as Shale: 
(Case numbers 89 and 117) 


Q = 0,10 to 0.20 
J = 0,15 to 0.20 
J' = 0,05 to 0.067 


Side friction of 25%. 

Using the average values of Q, J, and J' and side friction 
obtained for the various types of soil conditions, the variations of 
calculated ultimate resistances from load test results for the specific 
cases will be investigated. 

Hiley equation results are also shown with this comparison in 
Table CC, 

It should be noted that the values of ultimate resistance shown 
to be calculated by the Hiley formula were obtained by using values of 
C;, Co, and C3 recommended by Chellis (ref. 5) for use in the absence 
of specific information. No attempt was made to determine what changes 
in these coefficients would be necessary to achieve closer agreement 
with load test results. Comparison of Hiley formula ultimate re- 
sistance calculations has therefore been restricted to those case 
studies where complete data were not provided. Whether someone having 
experience with the Hiley formula could have justified use of different 
values for these coefficients and thereby achieved closer agreement 


with load tests is not known. 
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TABLE CC 


COMPARISON OF WAVE EQUATION SOLUTION 
WITH HILEY FORMULA AND LOAD TESTS 


Ultimate Resistance Per Cent Error 
Case Computed by Computed by Load Wave Hiley 
Number Wave Equation Hiley Formla Test uation Formula 
3 ho 52 5 ~11,1 -15.5 
3h 100 96 100 0.0 - 4.0 
Bo 103 92 100 3.0 - 8,0 
36 il 92 100 - 9.0 - 8.0 
67 60 25 60 0.0 -58.) 
a 6h, 80 6h 0.0 25.0 
(¢ 90 80 90 0.0 ~11.1 
76 50 60 50 0.0 20.0 
89 300 230 300 0.0 a23.3 
ge 300 10 300 0.0 -53.3 
pa3 100 122 100 0.0 2250 
Udy 122 90 17 4.3 -23.0 
15 133 86 182 ~26.9 =52.7 
1h6 97 72 91 6.6 -20.8 


161 103 96 167 -38,3 2.5 
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It must be remembered that the number of cases investigated 
here is small, and the values of point and side damping, ground quake, 
and side friction reported as giving good correlation may be replaced 
by more accurate values as more driving records and load tests are 
analyzed and as more basic information about the dynamic properties 
af soils becomes known. Furthermore, side resistance was taken as 
rectangular, or uniform distribution, along the side of the piles 
with about the top 10 feet of embedment discounted in each case, 
Different side distribution patterns would, of course, alter the values 
found above somewhat. 

On the other hand, the values reported are within the ranges 
that one might anticipate, and they do give good correlation as well 
as providing results generally more accurate than predictions from 


the Hiley formula, 
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Chapter VII 


RECOMMENDATIONS FOR FURTHER RESEARCH 


No attempt will, be made here to make an exhaustive and detailed 
list of further research possibilities, but several will be discussed 
briefly. Probably as many have already occurred to the reader, and 
certainly many more will present themselves to anyone beginning work 
in this general area. The main factors reqiiring further evaluation 
are ground quake, point and side damping, and distribution of side 
frictional resistance. Specifically, the recommendations ae as follows: 

A. Construct an apparatus in which model piles may be loaded 
both statically and dynamically. If economically feasible, the model 
pile may be instrumented and recording pressure cells embedded in the 
soil. Also it would be desirable to be able to erclose the entire 
soil-pile system so that a variable confining pressure could be 
applied. With this set-up the model pile could be encased to avoid 
side friction effects in the initial experiments. The model pile 
could then be struck by a small drop hammer and the set determined, 
This procedure could be repeated for different sized hammers and 
strokes, and the results compared with computer results for the test 
set-ups. The ultimate resistance could be determined from test load- 
ing statically an identical pile in the same controlled soil. In this 
way values of quake and damping might be bracketed more closely for 
various soilepile systems. 

B, Obtaining a picture of the interaction between pile and 


soil may prove useful. It is thought that this may be done for both 
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the static and dynamic cases. A photosensitive material could be 

used for a model pile which would be embedded in another photosensitive 
material representing the ground, The stress waves propagated by 
striking the model pile at its head by a small drop hammer could be 
captured in successive time intervals by using a high-speed camera. 
Even though such a study may give only qualitative results, it may 
assist in understanding how much of the ground is acting with a pile 
and the order of magnitude of stresses within this zone. Such informa- 
tion should prove helpful in evaluating both ground quake and the 
distribution of side frictional resistance. 

C, Study the problem from both an analytical and statistical 
viewpoint. Seek the cooperation of major construction corporations, 
government bodies, and others in obtaining as many good load tests 
as possible. For example, it is understood that the University of 
Michigan or the Bureau of Public Roads is on the verge of publishing 
anumber of very carefully controlled tests. Run a wide range of 
curves for various pile and hammer combinations as was done herein for 
the 53-pound bearing pile, and see if dynamic soil properties cainot 
be deduced from them which will give good correlation. If sufficient 
driving records and load tests become available, it may be advisable 
to write a computer program which will assist in matching up wave 
equation results with field data, 

D. A wide range of wave equation solutions may be run off 
for the common soil sampling rigs with a view toward correlating the 
information with the soil data and with pile driving and testing 


records. In this connection it should be noted that the pipe section 
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of the sampler should be broken down into one-foot increments and 
the time intervals of the calculations proportionately reduced in 
order to obtain stable calculations. It should also be kept in mind 
that since lateral support of the long pipe section is lacking in 
this case, some energy will be lost in transverse vibrations. How 
seriously this will affect the results is not known. 

E, The computer program included herein may be modified 
so that it can also handle the case of a long hammer, which should 
be represented as several weights in the mathematical model, and the 
case where the ram strikes the pile other than at its top. A. E. L. 
Smith illustrates how the mathematical model for these cases may be 
represented (ref. 28). 

F, More detailed study should be made of pile tip velocities 
at the time of maximum tip force, as determined from the computer re- 
sults, for a wide range of pile types, lengths, and types of hammers 
in order to establish a more valid basis for relating results of rapid 
loading tests made in the laboratory with pile driving action. 

G, With respect to the computer program, investigation should 
be made into the factors that cause instability in the computations. 

In a few computer results for this work, instability was encountered 
unexpectedly. These instances involved light piles, driven with 
relatively heavy hammers, or low values of ultimate resistance relative 
to the normal capacity of the pile, Investigations should be made to 
determine additional criteria for use in making up the input data so 

as to avoid instability. 

As previously indicated, these are but a few of the possibilities 
for further work; however, it is hoped that these listed will be suffi- 


cient to stimulate the interest of potential investigators. 
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CONCLUSION 


The prediction of a pile's ultimate static bearing capacity 
from its dynamic behavior during driving has been an elusive goal 
of civil engineers for many years. Hundreds of ingenious empirical 
dynamic formulas have been devised which may give good results under 
particular conditions of driving equipment, pile, and soil combina- 
tions; but the difficulty has been in selecting the proper formula 
to obtain reliable results. Furthermore, these dynamic formulas 
have limited ranges of applicability since each ignores important 
aspects of the problem. Because of these limitations, their validity 
as a general method of approaching the problem at hand is open to 
question. Civil engineers have recognized these shortcomings and 
have tended to treat the dynamic formulas as crude guides to the 
ultimate bearing capacity of a pile in granular soils, recognizing 
their complete inadequacy for piles driven in cohesive soils. 

The application of the wave equation to the pile driving 
problem, as envisioned here, represents an attempt to obtain a 
general method of analysis which will assist in predicting the 
ultimate bearing capacity of a pile from its observed behavior under 
the last hammer blow during driving subject to the limitations 
enumerated in paragraph E below. The wave equation method of analysis 
is also useful for studying stresses occurring during driving and 
for selecting appropriate pile driving equipment to meet economically 
specific field conditions. The major contribution by A. E. L, Smith, 


a pioneer in this work, was in devising a mathematical representation 
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of the pile driving problem suitable for digital computer solution, 
He then developed such a computer program, and his published work 
which followed emphasized the stress determination problem. He did, 
however, suggest that the method could be used in relating the static 
tearing capacity of a pile to its dynamic behavior, and he offered 
values of damping and ground quake which he believed to be accurate 
enough to make use of the method practical. 

To examine the possibility of correlating pile driving 
records with load tests by this method of analysis has been the 
principal objective of this thesis. Before making detailed case 
studies, it was necessary to examine the validity of the method 
itself, to review the soil engineering aspects of the problem, and 
to develop the necessary computer programs. These findings, together 
with the results from 2) case studies, are as follows: 

A, Validity of the mathematical model used. The model de- 
vised by A. E. L. Smith has been accepted as correct by those com- 
menting on his work, but it had not been proven to be equivalent 
to the wave equation. The formal proof of this equivalency has been 
provided in Chapter II. 

B. Resistance to @riving. It is shown in Chapter V that 
the expression used for the resistance of a pile to driving corresponds 
approximately to published experimental findings of others on the re- 
sistance developed in a soil due to high rates of loading. 

C. Computer programs, Although A E,. L. Smith developed 
computer programs for the wave equation method of analysis, they were 


not available as he had not published them. Computer programs are 
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published herewith which were developed independently, Their validity 
1s shown by comparing with both a manual solution for a sanple problem 
which utilizes all parts of the program and with results published by 
A, E. L. Smith. Complete instructions for usirg these programs are 
provided in Appendices A, B, and C so that someone unfamiliar with 
computer work may utilize them. 

The computer program developed for the Hiley formula is also 
subject to the limitations inherent in an empirical dynamic pile 
driving formula. 

D. Ground damping, quake, and side friction. From published 
experimental work by others, approximate values of damping to be used 
with the wave equation method of analysis for some soils have been 
derived in Chapter V. Starting with a range of quake reported by 
others, a range of values of damping were found which caused the 
solution of the wave equation method to correlate with load tests. 
This range of values of damping was compatible with the values derived 
from the published experimental work. Using the values of quake and 
damping thus established, values of uniformly distributed side fric- 
tional resistance were determined for various soil and pile combinations. 

E, Usefulness of the wave equation method for predicting 
static bearing capacity. Attempts to correlate 2) pile driving records 
with load tests have shown surprising results. It appears that the 
wave equation method of pile driving analysis may become an accurate 
method which may be used together with other factors in predicting the 


ultimate static bearing capacity of a pile, as indicated by the results 
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summarized in figure (657) and table (CC). Among these other factors 
which must be evaluated are the relaxation or set-up of the soil with 
time after driving, the group effect of piles, development of negative 
friction, consolidation of a clay layer beneath a pile tip, long-term 
changes in water table, corrosion or deterioration of the pile from 
other causes, and any other special problems which may be associated 
with a particular construction project. 

Furthermore, the problem has been treated as one-dimensional; 
i.e., lateral vibrations have not been included in the analysis. This 
should not, however, be a serious limitation on the method for most 
pile driving situations since a large proportion of the pile will be 
embedded as final penetration isreached. Hysteresis losses in the 
pile have not been included in the analysis, but their neglect is 
considered justified at this stage of development of the wave equation 
method, Since good correlation was possible using values derived 
from the published dynamic soil tests, it appears that the wave equation 
method is well suited for analysis of any size or type of pile being 
driven by any hammer in not only granular soils, but also cohesive soils, 

F, The dynamic interaction between pile and soil is incom- 
pletely understood and offers a fruitful area for further research. Even 
a qualitative picture of stress transfer should prove helpful. Addi- 
tional experimental work to isolate the parameters affecting ground 
quake, point damping, side damping, and distribution of point and side 
resistance to driving is needed as well as the evaluation of these 


factors. 


- 















tse ee nice pe Hy wa ak 
imu ~ — 6A eon Gl me Getta 
Chee oo ol | a cae 4 
a oe 4 
Te et om OI ew eee ERs ae AE 
Oden d ee oy 
SR ee Lene i ae Oe 

Se ae) er a) ie Ged ee 


mn cee ay) “eee | cee cee 78 aeitee em e658q) ée 


en MM) ae oo poe es ime eed pero te 





Sh Ee Te tle — ee ee Cle cee Og Or 
“oe Ae | ie of piety Cans) aut 
ieee anes 8 ee ee ee ee | 

eae Gir rr tA, » ~wede allie de oe ore i 
ihe ee eee « Mes —_— = & 


— eo ~~ i _2= "Vl ') (Vee tie aril een ak 


om sue- .* > fe © coms ob eerie oe . 
ae, «<6  » ite i (= |. oa cor Ap 


wn 2) La ee Pele! woe a. 


Sante ——. a) aes fm od ee “= —e if 


—_— 7 i. «= ers omit te Gael, ahem Or 





ee ie eto ce » -’ oe 6 md FIs, + Se ~~ ry ote ‘<-> 





ie (Gio ; senile. ol) tae a a) « wre * — '« et ae 





cme! ewe ai = ee 60 Ly OF Geta o! @ pore ie 


‘qar< 


Lue 


G. Correlation with many driving records and load tests 
should be made before taking the values reported herein as correct. 
Future investigators may use freely the computer programs and plotted 


results of calculations already made, and presented here, to facilitate 


such work. 
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APPENDIX A 
DETAILS OF THE "BASIC COMPUTER PROGRAM" 
WHICH USES THE WAVE EQUATION 

General: This program is written in Fortran language for the 
IBM 7090 computer, but it can be used equally well with the IBM 709 or 
IBM 709) computer. It has as its objective the detailed analysis of 
stresses, velocities, accelerations, displacements, and permanent set 
of the pile as these factors vary with time. For details of the theory 
and the mathematical model, refer to chapter III. 

Termination: Analysis is terminated automatically when the 
velocities of all the blocks of the mathematical model are simul- 
taneously negative or equal to zero. Negative velocities are taken 
as being in the direction from the bottom of the pile toward its top. 
Termination is also scheduled when the permanent set of the pile be- 
comes constant. Effects beyond either of these times are considered 
to be of secondary importance for purposes of this study. Generally 
these conditions will be met before calculations have been made for 
300 time intervals, and in most cases much sooner, The time interval 
is usually taken as 0.00025 or 0.00033 seconds because of the high 
velocity of stress wave travel in materials used for piles; therefore, 
this program investigates what happens during the first 0.1 second 
or less after the hammer hits the pile. In the usual case the pile 
Will achieve its permanent set well within this time, 

The program will also stop automatically if the velocity of 
either the pile cap or pile tip exceeds twice the velocity the ram 


had at the moment of impact. This feature is designed to help detect 
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instability (see Chapter III, C) in the calculations, but should not 
be relied upon entirely. The surest way of detecting instability is 
by plotting the displacements amd velocities of each block in the 
model against time and then examining the curves so obtained for very 
sharp peaks or discontinuities. Should such irregularities be present, 
the time interval may be reduced and the program rerun and rechecked 
for stability. The importance of this check cannot be over-emphasized. 
Should the program run the full 300 time increments, no 
indication is given in the output as to why the program stopped. In 
such a case the results should be examined especially carefully, for 
in all probability either instability was present or the program 
stopped before the final permanent set was reached. In such a case 
the variable DP may be redimensioned to a larger value (secund state- 
ment above statement number 207 in the beginning of the program), the 
statement immediately after statement number 500 should be changed to 
read DO 501 I =1,X where X will be the new dimension of variable 
DP, and the statement immediately after statement number 302 should 
be similarly changed to read DO 101 N = 3,X. In all other cases of 
program termination the reason for stopping will be printed out auto- 
matically in plain language; for example, "VELOCITY OF PILE CAP EXCEEDED 
TWICE THE RAM VELOCITY." 


Preparation of Input Data: All data is punched on standard 


data cards in the exact order indicated below and then placed immedi- 


ately after the * DATA card at the end of the source program. 
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CARD ORDER 


at 


The next M cards 


The next (M-1) cards 


The next M cards 


The next card 


The next card 
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DESCRIPTION 


Case number--any number up to 10 places may be 
assigned, and it should be punched within the 
first 10 places on the card, Blanks will be read 
as zeros; use no decimal points. 


If there is another case with more data following 
this set of data, punch +l in spaces 9 and 10 on 
this card; if this is the last case for this run, 
punch -] in places 9 and 10. 


In the first places of this card punch the num- 
ber of weights in the mathematical model. Within 
spaces 5 through 8 punch the number of weights 

in the mathematical model minus 1. Use no decimals. 


Punch within the first 10 spaces the decimal form 
of the time interval to be used in seconds, 


M is the number of weights in the mathematical 
model. One card is to be used for each weight. 
Punch the value in pounds of each weight in the 
first 10 spaces of each card. The first card is 
for weight 1 (the ram) and then consecutively 
down to the pile tip. Use a decimal point and 
the number can be located anywhere within the 
first 10 spaces. 


Punch within the first 10 spaces of each of the 
cards the spring constants beginning with the cap- 
block and extending to the top of the pile tip 
section of the model, The values are to be in 
pounds per inch. Use a decimal point and the num- 
ber can be located anywhere within the first 10 
spaces, 


Punch within the first 10 spaces of each card the 
soil spring coefficients beginning with block Il, 
At least the first 2 or 3 will be zero normally. 
Use a decimal point and the number in pounds per 
inch can be placed anywhere within the first 10 
spaces. 


Punch with a decimai point anywhere within the 
first 10 spaces on this card the value of the 
capblock coefficient of restitution, 


Punch with a decimal point anywhere within the 
first 10 spaces on this card the value of the 
coefficient of restitution of the head packing 


ie © 4m < 
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in the case of a concrete pile. For all other 
cases without head packing punch 1,00, 


The next card Punch with a decimal point anywhere within the 
first 10 spaces the velocity of the ram at the 
instant of impact in feet per second. 


ine next card Punch with a decimal point anywhere within the 
first 10 spaces the value of ground quake in inches. 


The next card Punch with a decimal point anywhere within the 
first 10 spaces the value of point damping to be 
used. 

The next card Punch with a decinal point anywhere within the 
first 10 spaces the value of side damping to be 
used, 

The next card If the pile cap is allowed to transmit tension, 


enter a +1 in spaces 9 and 10; if not, enter a 
-l there. Use no decimals. 


The next card if side friction is present, enter a +1 in spaces 
9 and 10; if not, enter a -1 there. Use no 
decimals, 


Running the Program: The data deck of cards just described 
1s placed immediately after the source program deck and is ready to be 
turned over to a machine operator for processing. 

Results: The computer will clearly write out in the output 
the input datas the number of each time interval; and identify the dis- 
placements, velocities, permanent sets, forces, and ground resistances 
acting on each block. The results are printed in decimal form with a 
multiplier notation, For example, 0.2207 E06 means 0.2207 x 10° or 
220,700. 

The Program: The complete program is presented so that any 
who wish to use it may do so freely. Each line of the program is punched 
on a separate card and must be assembled in the exact order shown. Any- 


one preparing data cards from this program should verify painstakingly 





oe be Se Oe oe ey ee 
oo.) Ott oe Gere eee 16 gee 


—- Ge hee te ro fe 
a" oe ee’ @ ees! 
: een “ep '  @ @! eee) «gp 





—* «'S _ oe ee) \ * ee 
—s i al i om; 7 —P -& =, *f ter 


WJ, eee Pre’ wer «a _. 
Ss te. “ae ae 3 anes ¥ ial 





ms} « —— * = es 5 “Hoe 
' e-ae > oual't — aha 4 jiewi 
——_™ 


—————— —_———— *: aus e @) ae a" ome Aj 


a? + -_-teern © | he 


ee —-=_ tH te 





~~ -* ae eee e) cee OS C. a 


oom! fs oti ~ e ot a feet 







Ah Fe ah le: ae See 
5 wv : - Tm ule lh tl em] 





‘< = =) aa ) «2 TT 


ae. The os ie ———- wi —J ae, 
- . ° ‘ eA Sa 4 a ay ae 
S vo ° .s te mg) x ¢ mtg © . 
—— 5 2 i iL mm." 
: - o Abad 
' —_/ —— — 


ize 


that no mistakes have crept in during punching, for the slightest dis- 
crepancy can lead to totally erroneous results. For that reason those 
preparing cards from the program given here are cautioned to try it 
out with test cases before relying on the results. 

The notation used in the program is the same as was used in the 
derivation of the basic equations except as dictated by Fortran require- 
ments. The more important variations in terminology are: 

S instead of K for spring constants 

SP instead of K! for soil spring coefficients 


RES] instead of ey for coefficient of restitution for the 


capblock 

RES2 instead of en for coefficient of restitution for 
head packing 

Z instead of J for point damping factor 


ZP instead of J! for side damping factor 
Additionally, the distinction between the present time interval, the 
last time interval, and the second previous time interval is made by 
use of the letters L, M, and S respectively. For example, in the case 
of velocities, 

VL is used instead of V 

VM is used instead of v 

VS is used instead of v* 
Other variables introduced in the program are defined as they occur 


within the program or have been explained under preparation of data 


above. 
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SID REESE C2165 CIVIL ENGR GWAVE EQ 5 SO 
% XEQ 
or LABEL 


s WAVE EQUATION 
DIMENSION DS(100), DM(100), DL(100),VS(100), VM(100), VL(100) 
DIMENSION DPSM(100), DPSL(100) 
DIMENSION S(100), W(100), DP(300), R(100), SP(100),CS(100) 
DIMENSION CM(100), CL(100), FS(100), FM(100), FL(100), DPSS(100) 
207 DO 500 I = 1,100 
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FS(I) 
FM(T) 
FL(I) 
500 DPSS(I O 
DO 501 a - 300 
E@in DE(L) = 
READ eis Sar: 5,27,NCASE, ITURN 
27 FORMAT (110/110) 
O READ INPUT TAPE 5, 70,M,1IM1,T, (W(L),I=1,M), (S(I),I=1,00) 
L , (SPQ), 1=1,™), RESL, RES2, v One. 2P, CAP, RUB 
70 FORMAT (211 /(F10. no) 
WRITE OUTPUT TAPE 6,30,NCASE 
30 FORMAT (12H CASE NUMBER I5) 
WRITE OUTPUT TAPE 6,28,(I,W(I),S(I), s(t), T=1,M) 
28 FORMAT (SH DATA//58H Me S 
1 SP/(Ih, 3F20.2)) 
WRITE OUTPUT TAPE 6 ,29,RESL,RES2,ZP,2,Q,V 
29 FORMAT (36H COEF. OF REST, OF CAPBLOCK (RES1) = Fh .2/ 
1 36H COEF, OF REST. OF PILECAP (RES2) = Fh,2/ 
2 36H SIDE DAMPING FACTOR (ZP) F.2/ 
3 36H POINT DAMPING FACTOR (2) Bias 
)} 36H GROUND QUAKE (Q) = Fh.2/ 
5 36H INITIAL RAM VELOCITY (V) F6.2) 
WRITE OUTPUT TAPE 6,31,T,CAP, RUB 
31 FORMAT (36H TIME INTERVAL = F8.6/ 
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36H CAP = F6,2/ 

36H RUB = F6.2) 

WRITE OUTPUT TAPE 6,1 

FORMAT (8H M DL VL DP 
FL R//) 

DS(1L) = V#12.«2 

VS(1) = V +(-DS(1)#S(1) )#T%32.2/(1) 

VS(2) = (DS(1)#S(1) )¥T+32,2/N(2) 

DM(1) = DS(1) + VS(1)+12.s7 

DM(2) = VS(2)+12,%T 

CM(2) = DM(2) 

VM(1) = VS(1) -((DM(1)-DM(2 ) )#S(1) )+T%32.2/W(1) 

VM(2) = VS(2)+((DM(1)-DM(2) )«S(1)-(DM(2)-DM(3))#S(2) )#T#32,2/W(2) 

vM(3) = VS(3) +((DM(2)-DM(3))*S(2)-R3 )*T#32,2/w(3) 


CM(1) = DM(1) - DM(2) 

CHECK = 2,.% V 

WRITE OUTPUT TAPE 6,301,DS(1), VS(1) 
FORMAT (SH N= 1,68 1 2E15.4//) 
WRITE OUTPUT TAPE 6, 302,DM(1),VM(1),DM(2), VM(2) 
FORMAT (QHN = 2,68 1 2E15.4/14H 2 2E15.4//) 
DO 101 N = 3, 300 

WRITE OUTPUT TAPE 6, 303,N 

FORMAT (5H = 13) 

MLESS1 = M-1 

DO 130 I = 1,M 

DL(I) = DM(I) + VM(I)#12.%T 

DO 131 I = 1,MLESS1 

CB(L) = DL(I) = DL({i+tl) 

MCDL(M)) 12,12,14 

DE = 0.0 

GO TO 2h 

IF(DL(M) - Q) 16,16,18 

DE = 0.0 

GO TO 2h 

VALUE = DL(M) - Q 

IF (VALUE-DE) 22,22,20 

DE = VALUE 

GO TO 2h 

DE = DE 

DP(N) = DE 

R(M) = (DL(M)-DP(N))*SP(M)#(1.+Z*VM(M) ) 
IF (DP(N)) 25,25, 26 

IF (DP(N) - DP(N-1)) 197,197, 25 
CONTINUE 

VAL = S(1)#CL(1) 

IF(CL(1)-CM(1)) 5,5,h 

Pa) = VAL 

GO TO 33 

CL(1) = CM(1) 

FL(1)=VAL/(RESL##2)=(1, /(RESH4+2 )=1, )##2(1)*CL(1) 
VAL2 = S(2)#CL(2) 
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FL(2) = VAL 

GO TO 37 

CL(2) = CM(2) 

FL(2) = VAL2/(RES2#%2)~(1, /(RES2##2 )-1, )#S(2)#CL(2) 
IF(CAP) 38, 36,36 

on. 39,36, 36 

mie) = 

ae 

PERPL(T)) 8,133,133 

Pi) = 0,0 

DO 132 I = 3,MLESS1 

FL(I) = CL(I)xS(I) 

IF (RUB )h9,09, bh 

DO 48 I = 3, MLESS1 

DPSL(I) = DESL) 

CHANGE = aD) e 

SUM + DL(I) + Q 

IF (DPSL(I)-CHANGE )h5, 4.6, 6 

DPSL(I) = CHANGE 

=: 48,47 

DPSL(I) = 

R(I) = iL) - DPSL(I) )*SP(I)*(1.+ZP*VM(I)) 
DO 50 I =1,M 

ving) = VM(L) + (PFL(I-1L) - FL(I) - R(I))#T%32.2/w(I) 
CONTINUE 

DO 100 I =1,M 

WRITE OUTPUT TAPE 6 , 304, 1 Ta VL(I),DPSL(I),FL(I),R(1I) 
FORMAT (11H 13,5E1 Ly) 

CONTI NUE 

DO 75 K=1,M 

ee, 

dct (VL(K)) 75,75, 800 

CONTINUE 

GO TO 199 

CONTINUE 

IF(VL(2) - CHECK)99,99,190 
IF(VL(M)-CHECK )102 ,102,19) 

DO 98 K=1,M 


DS(K) = DM(K) 
DM(K) = DL(K) 
DL(K) = 0,0 
VS(K) = VM(K) 
VM(K) = VL(K) 
VL(K) = 0,0 
CS(K) = CM(K) 
CM(K) = CL(K) 
CL(K) = 0.0 
FS(K) = FM(K) 
FM(K) = FL(K) 
FL(K) = 0.0 
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DPSS(K) = DPSM(K) 

DPSM(K) = DPSL(K) 

DPSL(K) = 0.0 

WRITE OUTPUT TAPE 6,305,DP(N) 
FORMAT (1hH E15.) 
CONTINUE 

GO TO 196 

WRITE OUTPUT TAPE 6,191,N 


hou W 


191 O FORMAT (6H VELOCITY OF PILE CAP EXCEEDED TWICE THE RAM VELOCITY W 


d. 


19k; 
195 0 
a 


a 
198 


7? 
200 
196 
205 


+ 


HEN N WAS 13) 

GO TO 196 

WRITE OUTPUT TAPE 6,195,N 

FORMAT (6H VELOCITY OF PILE TIP EXCEEDED TWICE THE RAM VELOCITY W 
HEN N WAS I3) 

GO TO 196 

WRITE OUTPUT TAPE 6,198,N 

FORMAT (26H DP BECAME CONSTANT AT N = 13) 

GO TO 196 

WRITE OUTPUT TAPE 6,200,N 

FORMAT (52H ALL VL WERE SIMULTANEOUSLY NEGATIVE OR ZERO AT N = I) 
IF (ITURN )205 , 205 , 207 

CALL EXIT 

END 

DATA 
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APPENDIX B 


DETAILS OF THE "VARY Ru PROGRAM" WHICH USES THE WAVE EQUATION 


General. This program is identical to the BASIC PROGRAM except 
that it provides the capability of varying certain input data in order 
to obtain wave equation solutions for a range of input variables. The 
program provides for a maximum of 500 time cycles, instead of 300, to 
allow for the use of smaller time increments and a greater number of 
pile segments. Otherwise, the termination features are the same as for 
the BASIC PROGRAM. 

Program Features. The program is designed with the following 
features that are different from or in addition to those of the basic 
program, 

1. The ultimate resistance (Ru) is varied over a preselected 
range of values in any increment desired. The set is computed for 
each value of Ru. The initial value of Ru, the increment of increase 
of Ru and the number of times the increment is applied are provided as 
input data. 

2. the effective length of embedment of the pile can be varied 
and is determined by input data. This is the length of the pile over 
which the resistance to driving acts. 

3. Five patterns of resistance distribution between point bear- 
ing and side friction can be applied for computation. These patterns 
vary from full point bearing to full side friction with intermediate 
distributions of 75, 50, and 25 per cent point bearing with the balance 


distributed as side friction. 
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4, The method of distributing side friction can be rectangu- 
lar or triangular as desired. Input data determine the selection of 
the method. 

5. The capability of recycling the entire program up to 
three times with the new values of any variable is provided. This 
feature was used to provide new values of Q and J, although other 
variables could be used. This capability is provided with acontrol 
Statement, and the new values of the variables to be changed must be 
put into the program itself, a somewhat awkward arrangement. 

6. Any number of sets of data can be used in a single run, 
as provided for in the input data. This is also a feature of the 
basic program. 

Input Data, The input data is the same as for the basic 
program except for the additional data required. 

CARD ORDER DESCRIPTION 

2a After the second card of the basic program, 

areas of the pile head, center and tip in square 
inches are punched on one card in the first 
three ten-column fields, Values are punched with 
a decimal point. 

After the last card In the first four columns punch 1 to indicate 

of the basic program triangular distribution of resistance on the side 
or -1 to indicate rectangular distmbution. 

Next card In the first 10 columns punch the initial value 
of Ru in pounds, and in the second 10-column 
field the increment of Ru in pounds. A decimal 
point must be used, These are variables RU1 
and ADD. 

Next card In the first four columns punch the number rep- 
resenting the number of times the increment of 


Ru is to be applied; this establishes the maximum 
value of Ru. In the second four columns (5 - 8) 





Next card 


Next and last card 


Output Data. 


ing output is provided: 


le 


punch the number of blocks in the mathematical 
model to which ground resistance is not to be 
applied; this establishes the effective embedded 
length of the pile. These variables are ICYCLE 
and NOLOAD and are punched without a sign or 
decimal point, 


A number of 1 through 5 is punched in each of the 
first three four-colum fields for the variables 
LA, LB, and LC respectively. These numbers pro- 
vide the values for the indexing parameters of 

a DO statement that determines the distribution 
of resistance between point and side. 


A number with a value between 1 through ); is 
punched in the first four columns. This number 
controls the number of times the program is to 

be repeated with new variables. No sign or deci- 
mal is used. This variable is called IRES. If 
the number 1 is used the program will cycle the 
full three times, and if is used it will stop 
after all computations using the initial data 

are made, 


In addition to listing all input data, the follow- 


1. Value of Ru in pounds for the particular computation. 


2. oet in blows per inch, 


3. Compression and tension stresses at the head, center and 


tip of the pile, in pounds per square inch, 


4. The number of time cycles run, 


Use of the Program, This program is designed to provide the 


data necessary to plot Ultimate Resistance vs, Set curves for various 


values of Q, J,and J', and for one or more of the five different dis- 


tributions of resistance between point and side. The values of Q, J, 


and J' are initially provided in the input data. If additional values 


are desired for computations, they are provided through the use of the 


"computed-go-to" statement following statement number 600, The 
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particular values must be placed in the program following statements 
700, 7Ol,and 702. This feature could easily be expanded if desired. 

If it is desired to use all five distributions of resistance between 
point and side, the values of LA, LB, and LC must be 1, 5, and 1 re- 
spectively. When the index variable of the DO statement (five state- 
ments after 777) takes different values, the distribution of resistance 


is aS indicated below: 


Value Distribution 
1 All point bearing 
2 75 per cent point bearing, 25 per cent friction . © 
3 50 per cent point bearing, 50 per cent friction 
hh 25 per cent point bearing, 75 per cent friction 
5 All side friction 


By assigning proper values to LA, LB, and LC, the desired distributions 
are obtained in accordance with the operation of a DO statement as ex- 


plained in reference (21). 
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—_ C2165 CIVIL ENGR G WAVE EQ 15 250 VARY RU 

Q 

LABEL 

WAVE EQUATION VARY RU 101 

DIMENSION DS(100), DM(100), DL(100),VS(100), VM(100), VL(100) 

DIMENSION DPSM(100), DPSL(100 ) 

DIMENSLON S(100) w(100), DP(500), R(100), SP(100),CS(100) 

DIMENSION CM(100), CL(100), FS(100), FM(100), FL(100), DPSS(100) 

READ INPUT TAPE 5,27,NCASE, ITURN 

FORMAT (110/110) 

READ INPUT TAPE 5,71, AREAL, AREA2, AREA3 

FORMAT ( 3F10, 1) 

READ INPUT TAPE 5, 70, M,IML,T, (W(I),I=1,M),(S(I),I=1,1M) 
»(SP(I), 1=1,M), RESL, RES2, y, QZ. ve CAP, RUB 

FORMAT (21h/(F10.h)) 

READ INPUT TAPE 5,60), MUD 

FORMAT(I) ) 

READ INPUT TAPE 5,72,RU1,ADD,ICYCLE,NOLOAD 

FORMAT (2F10, 1 /2Th ) 

READ INPUT TAPE 5,73,LA,LB,LC,IRES 

FORMAT (314/Ih) 

WRITE OUTPUT TAPE 6,30,NCASE 

FORMAT (12H CASE NUMBER I5 ) 

WRITE OUTPUT TAPE 6,28, (I,W(I),S(I), ihe T=1,M) 

FORMAT (5H pata//S8H M S 

SP/(Ih, 3F20.2)) 

WRITE OUTPUT TAPE 6,29,RES1,RES2,ZP,2,Q,V 

FORMAT (36H COEF. OF REST. OF CAPBLOCK (RES1) = F.2/ 

36H COEF, OF REST, OF PILECAP (RES2) = Fh.2/ 

36H SIDE DAMPING FACTCR (ZP) Fh.2/ 

36H POINT DAMPING FACTOR (Z) Fu.2/ 

36H GROUND QUAKE (Q) = Fh.2/ 

36H INITIAL RAM VELOCITY (V) F6.2) 

WRITE OUTPUT TAPE 6,605,MUD 

FORMAT (6H MUD = 16) 

WRITE OUTPUT TAPE 6,31,T,CAP,RUB 

FORMAT (36H TIME INTERVAL = F8.6/ 

36H CAP = F6.2/ 

36H RUB = F6.2//) 

WRITE OUTPUT TAPE 6,608,RU1, ADD, ICYCLE, NOLOAD 

FORMAT (5H RU = E20. , 10H ADD = E20.4/9H ICYCLE = 16, 

13H NOLOAD = I6) 

WRITE OUTPUT TAPE 6,609,LA,LB,LC,TRES 

FORMAT (5H LA = Ih, 10H LB = oun 10H LC = Ih,10H TRES = 

Th) 

CHECK = 2,.% V 

CONTINUE 

RU = RU1 

DO 600 IAM = 1,ICYCLE 
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ITEM = M - NOLOAD 
AITEM = ITEM 
DO 602 J = LA,LB,IC 
FMAXL 
FMAX2 
FMAX3 
FTENI 
FTEN2 
PTEN3 
DO 500 I = 1,100 
DS(I) 
DM(I) 
DL(I) 
VS(I) 
VM(I) 
) 
( 
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Re 
Vx12 .*T 
V +(-DS(1)*S(1L) )xTx32.2AN(1) 
(DS(1)#*S(1))#Tx32,2/u(2) 
DS(L) + VS(122.xT 
VS(2 )612 xT 
DM(2) 
VS(1) -((DM(1)-DM(2) )#S(1) )T#32.2/W(1) 
VS(2)+( (DM(1)-DM(2 ) )S(1)-(DM(2)-DM( 3) )*S(2) )#T#32 .2/w(2) 
VS(3)+( (DM(2)-DM(3) )#S(2)-R(3) )#T#32.2/N(3) 
DM(1) = DM(2) 
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DO sor I = 1,500 
) 
) 
) 
) 
) 
) 
) 


Q 
22 
vr 
f= 
nee? Ne” 
non 


f= J 
PART = (5,-AJ)/\. 
SP(M) = PARTxRU/Q 
SIDE = RU - PART#RU 
LOAD = NOLOAD + 1 

IF (MUD )650, 651,652 

DO 601 N = LOAD, TM 
SP(N) = SIDE /(ALTEMxQ ) 
SP(M) = SP(M) + SIDE/(ALTEMxQ) 
GO TO 651 

DO 603 N = 3,1M 

NT = N-2 
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fel = Q2xNT . 1 


IM2 = M2 
ANTL = NTL 
Aivi2 = IM2 


SP(N) = SIDE*ANT1/( (AIM2+*2 )xQ) 
SP(M) = SP(M) + SIDEx(2.sAIM2-1. )/( (ATM2s2 )#Q) 
CONTINUE 

DO 101 N = 3,500 

MLESS1 = 

DO 130 I = 

DL(I) = a + VM(L)x12."9 

DO 131 I = 1,MLESS1 

CL(I) = DL(Z) - DL(I#+1) 
IF(DL(M)) 12,12,14 

DE = 0,0 

GO TO 2h 

IF(DL(M) - Q) 16,16,18 

DE = 0.0 

GO TO 2h 

VALUE = DL(M) - Q 

IF (VALUE-DE) 22,22,20 

DE = VALUE 

GO TO 2h 

DE = DE 

DP(N) = DE 

R(M) = ce “SP(M)#(1.+Z*VM(M) ) 
TF (DP(N)) 25,25,26 

IF (DP(N) - DP(N-1)) 197,197, 25 
CONTINUE 

VAL = S(1)#CL(1) 
TF(CL(1)-CM(1)) 5,5,4 

FL(1) = VAL 

GO TO 33 

Cee] CM) 

FL(1 )=VAL/(RESLiv¢2 )-(1. /(RESL4¢2 )-1. )*S(1)#CL(1) 
VAL2 = S(2)*CL(2) 
IF(CL(2)-CM(2)) 35,35,34 


FL(2) = VAL2 
GO TO 37 
CL(2) = CM(2) 


FL(2) = VAL2/(RES2*%2 )-(1./(RES2%"2 )-1, )#5(2)#CL(2) 
IF(CAP) 38,36,36 

> 39, 36, 36 

FL(2) = 

pl 

a) =? 
FL(1) = 

DO 132 ae . MLESS1 

FL(I) = oh “S(L) 

IF (FL(2) )213,211, 210 

IF (PIEN1-FL(2))212, 211,21) 
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FTENL = FL(2) 

GourTO 211 

IF (FMAX1-FL(2))215,211, 211 
FMAX1 = FL(2) 

CONTINUE 

MC = M/2 +1 

IF(FL(MC ) )216, 217, 218 

IF (FTEN2-FL(MC))217, 217,219 
FTEN2 = FL(MC) 

COMTO 217 

IF (FMAX2-FL(MG ) )220, 217,217 
FMAX2 = FL(MC) 

CONTINUE 

IF (FL(IM ) )221, 222, 22h 

IF (FTEN3-FL(IML) )222, 222, 223 
FTEN3 = FL(IM) 

GO TO 222 

IF (FMAX3-FL(IML) )225, 222,222 
FMAX3 = FL(IML) 

CONTINUE 

TF(RUB)L9,H9, 4h 

DO h8 I = 3, MLESSI 

DPSL(I) = DPSM(I) 

CHANGE = DL(I) - Q 

SUM = DL(I) + Q 

IF (DPSL(I )-CHANGE )5 , 6,6 
DPSL(I) = CHANGE 

ee uN6, 48,47 


DPSL(I) = 
R(I) = (BU(E) - DPSL(I))#SP(I )(1. +ZP*VM(I ) ) 
DO 501 =1,M 


VL(I) = vM(Z) + (FL(I-1) - FL(I) - R(T))#032.2/W(1) 
CONTI NUE 

DO 75 K = 1,M 

X=, 

IF (VL(K)) (ay >, 000 
CONTINUE 

GO TO 199 

CONTINUE 

IF(VL(2) - CHECK)99,99,190 
IF (VL(M)-CHECK )102, 102 "19h 
DO 98 K = 1,™ 

DS(K) = DM(K) 


DM(K) = DL(K) 
DL(K) = 0.0 
VS(K) = VM(K) 
VM(K) = VL(K) 
VL(K) = 0.0 
CS(K) = CM(K) 
CM(X) = CL(K) 
Gein) = 0,0 
FS(K) = FM(K) 
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Pew, = 0,0 
DPSS(K) = IPSM(K) 
DPSM(K) = DPSL(K) 
DPSL(K) = 0.0 
CONTINUE 

GO TO 196 


WRITE OUTPUT TAPE 6,191,N 


O FORMAT (64H VELOCITY OF PILE CAP EXCEEDED TWICE THE RAM VELOCITY W 
1 HEN N WAS I3) 


GO TO 196 
WRITE OUTPUT TAPE 6,195,N 


O FORMAT (6)H VELOCITY OF PILE TIP EXCEEDED TWICE THE RAM VELOCITY W 
1 HEN N WAS 113) 


0 


GO TO 196 

WRITE OUTPUT TAPE 6,198,N 

FORMAT (26H DP BECAME CONSTANT AT N = T3) 

GO TO 196 

WRITE OUTPUT TAPE 6,200,N 

FORMAT (52H ALL VL WERE SIMULTANEOUSLY NEGATIVE OR ZERO AT N = Ih) 
CONTINUE 

WRITE OUTPUT TAPE 6,607,RU 

FORMAT (27H ULTIMATE RESISTANCE (RU) = E15.) 
WRITE OUTPUT TAPE 6,330,Q,W(1),V 

FORMAT (4H Q = E10.2,7H Wi) = EIS. WH V = £15.) 
WRITE OUTPUT TAPE 6, 606, (SP(L), L=1 .M) 

FORMAT (17H NEW VALUES OF SP/(BE15.))) 

WRITE OUTPUT TAPE 6,305,DP(N) 

FORMAT (7H SET = E10.) 

BLOW = 1./DP(N) 

WRITE OUTPUT TAPE 6,805, BLOW 

FORMAT (66H 

BLOWS/IN= E10.3) 

FTEN] = FTENI/AREA1 

FTEN2 = FTEN2/AREA2 

FTEN3 = FTEN3/AREA3 

WRITE OUTPUT TAPE 6,307,FTEN1,FTEN2,FTEN3 
FORMAT(30H MAXIMUM TENSION AT HEAD = ELS .y/ 
30H MAXIMUM TENSION AT CENTER = E15.)/ 

30H MAXIMUM TENSION AT TIP = E15.) 

FMAX1 = FMAX1/AREAL 

FMAX2 = FMAX2/AREA2 

FMAX3 = FMAX3/AREA3 

WRITE OUTPUT TAPE 6,306, FMAX1, FMAX2 ,FMiX3 


it il 


FORMAT (26H MAXIMUM FORCE AT HEAD = £15.),/ 
26H MAXIMUM FORCE AT CENTER = E15.),/ 

26H MAXIMUM FORCE AT TIP = £15.//) 
CONTINUE 

RU = RU + ADD 


CONTINUE 
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GO TO (700,701, 702,703), IRES 
CONTINUE 

7a=mOn 20 

ZP= 0.067 

TIRES = 2 

GO TO 777 

CONTINUE 

Q = 0,20 

TRES = 3 

GOenO. 777 

CONTI NUE 

Z= 0,15 

ZP = 0.05 

TIRES = 

GO TO 777 

CONTINUE 

IF (ITURN )205 , 205 , 207 
CM EXIT 

END 

DATA 
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APPENDIX C 
DETAL LS OF THE "RESEARCHER COMPUTER PROGRAM" 
WHICH USES THE WAVE EQUATION 

General: This program is written in Fortran language for the IBM 
7090 computer, but it can be used equally well with the 709 or 709) 
computer, It has as its objective the systematic production of data 
fron which set vs. ultimate ground resistance curves can be plotted 
for any case. Its chief advantage lies in the fact that for every set 
of data run, the vrogram automatically begins with a value of ground 
resistance of 20 tons and finds the set associated with all the load 
carried by the pile tip, 75% of load carried by the pile tip, 50% of 
load carried by the pile tip, 25% of load carried by the pile tip, and 
finally all the load carried by friction along the sides of the pile. 
It then assumes an increase of ground resistance of 20 tons and repeats 
the calculations. It automatically assigns new values of ground re- 
sistance in this manner up to and including 280 tons. Essentially, 
then, it progresses through the basic program five times at each value 
of ultimate ground resistance, or 70 times for each set of data, It 
also writes out the maximum compression and tension stresses in pounds 
per square inch for each of these 70 conditions at the head, mid~length, 
and tip of the pile. The side friction distribution may be chosen as 
rectangular or triangular in shape. In both of these cases the top 10 
feet of the pile is considered free of frictional resistance. 

Termination: Termination of each of the 70 conditions for each 
set of data is programmed exactly like the basic program. The reason 


for each termination is written out, and the computer immediately moves 


aye 
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to the next set of conditions. No over-all termination between sets 
of data should occur. 

stability: Besides the checks to detect instability in the 
calculations as described under the basic program, a further check 
is obtained when plotting set vs. ultimate ground resistance. This 
plot should also produce smooth curves. If not, a set of data for 
each suspicious point can be prepared and run with the basic program 
to obtain a detailed analysis so that evaluation can be made of the 
validity of the results. leduced time intervals can be used to elin- 
inate the instability, or smaller pile lengths in the mathematical 
model together with reduced time intervals may be used. 

Preparation of Input Data: All data is punched on standard 
cards exactly as described in the basic program except that two addi- 
tional cards are required far this program. 

The first extra card is inserted between cards 2 and 3. I 
the first 10 spaces of this new card are punched, with a decimal point, 
the cross sectional area of the pile at its top, within spaces ll 
through 20 is punched in the same manner the pile cross-sectional area 
at mid-length of the pile, and similarly the area of the pile tip in 
spaces 21 through 30. 

The second extra card is included as the very last card of 
data. If side distribution is rectangular, punch -1 in spaces 3 and 
hs; if side distribution is to be triangular, punch +l in spaces 3 and kh. 

In all other respects, data is prepared exactly as described 


under instructions for the "basic computer program." 
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Running the program: The data deck of cards just described 
is placed immediately after the program deck (after the * DATA card) 
and is then ready to be processed by a machine operator. 

Results: The input data will be repeated in the output 
information. The time cycle that the calculations were terminated 
for the first condition and the reason therefor will first appear. 
The ultimate ground resistance; the soil spring coefficients; set; 
and the maximum tension and compression at the pile head, mid-length, 
and tip will be written out for each condition. The results are 
printed in decimal form with a multiplier notation. For example, 
0.1511 EOl means 0.1511 x 10° = 1.511. 

The Program: The complete program is presented so that any 
who wish to use it may do so freely. Each line of the program is 
punched on a separate card and must be assembled in the exact order 
given below. Extreme care must be used in punching to prevent errors 
from creeping in, and for that reason it is recommended that test 
cases be run to check the program before relying on its results. 

The notation used is as described under this section of 


instructions for the "basic computer program." 
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REESE C2165 CIVIL ENGR G WAVE RQ 4 100 RESEARCHER 
XEQ 

NOBIN 

LABEL 

WAVE EQUATION VARY RU 

DIMENSION DS(100), DM(100), DL(100),VS(100), VM(100), VL(100) 
DIMENSION DPSM(100), DPSL(100) 

DIMENSION S(100) W(100), DP(300), R(100), SP(100),CS(100) 
DIMENSION CM(100), CL(100),FS(100), FM(100), FL(100),DPSS(100) 
READ INPUT TAPE 5 27, NCASE, TTURN 

FORMAT (I10/I10) 

READ INPUT TAPE 5,71, AREA], AREA2, AREA3 

FORMAT(3F10.)) 


O READ INPUT TAPE 5, 70, M,IM,T, (W(I),I=1,M), (S(1),I=1, IM) 


1 


fi 


i 
2 
3 
n 
5 


1 
2 


, (SP(I), I=1,™), RESI, RES2, y, Q,2Z, ZP, CAP, RUB 
FORMAT (21 /(F10.h)) 
READ INPUT TAPE 5,60, MUD 
FORMAT(IbL ) 
WRITE OUTPUT TAPE 6,30,NCASE 
FORMAT (12H CASE NUMBER I5) 
WRITE OUTPUT TAPE 6,28, (I,W(I),S(Z), = I=1,M) 
FORMAT (5H pata//58H M . 
Bey (ilig3h20.2)) 
WRITE OUTPUT TAPE 6,29,RES1,RES2,ZP,Z,Q,V 
FORMAT (36H COEF, OF REST, OF CAPBLOCK (RES1) = Fh.2/ 
36H COEF. OF REST, OF PILECAP (RES2) = Fh.2/ 
36H SIDE DAMPING FACTOR (ZP) 
36H POINT DAMPING FACTOR (Z) 
36H GROUND QUAKE (Q) = Fh.2/ 
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36H INITIAL RAM VELOCITY (V) = F6.2) 
WRITE OUTPUT TAPE 6,605,MUD 
FOAMAT (6H MUD = I6) 
WRITE OUTPUT TAPE 6,31,T,CAP,RUB 
FORMAT (36H TIME INTERVAL = F8.6/ 
36H CAP = F6.2/ 
36H RUB = F6,2//) 
CHECK = 2.% V 
= },0000, 
DO 600 I = 1,14 
ITEM = M-3 
AITEM = ITEM 
DO 602 J = 1,5 
FMAX1 = 0.0 
FMAX2 = 0.0 
FMAX3 = 0.0 
FTEN1 = 0.0 
FTEN2 = 0.0 
FTEN3 = 0,0 


DO 500 I = 1,100 
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DS(I) = 0.0 

DM(I) = 0,0 

DL(I) = 0.0 

Va(I) = 0.0 

VM(I) = 0.0 

VL(I) = 0.0 

DPSM(I) = 0.0 

DPSL(I) = 0,0 

Ren = 0.C 

SP(I) = 0.0 

CS(I) = 0.0 

CM(I) = 0.0 

CL(I) = 0.0 

bower) = 0.0 

FM(I) = 0.0 

fel) = 0,0 

DPSS(I) = 0,0 

DO 501 I = 1,300 

He) = 0.0 

DS(1) = V¥12.+7 

VS(1) = V +(-DS(1)#S(1) )xT#32.2/w(1) 

VS(2) = (DS(1)#«S(1))#Tx32.2/W(2) 

DM(1) = DS(1) + VS(1)#12.#T 

DM(2) = VS(2)s12 22 

CM(2) = DM(2) 

VM(1) = VS(1) -((DM(1)-DM(2) )2S(1))#T#32.2/W(1) 
VM(2) = VS(2)+((DM(1)=-DM(2 ) )%*S(1)-(DM(2 )-DM( 3) )S(2) )xT#32.2/W(2) 
VM(3) = VS(3)+((DM(2)-DM(3) )*S(2)-R(3) )Tx32.2/W(3) 
CM(1) = DM(1) - DM(2) 


PART = (5,-Ad)/k. 

SP(M) = PART+RU/Q 

SIDE = RU - PART#RU 

IF (MUD )650, 651,652 

DO 601 N = “h 10 

SP(N) = SIDE/( ALTEM*Q) 
SP(M) = SP(M) + SIDE/(ATTEMQ) 


GO TO 651 

DO 603 N = 3,1 
Ni anl=2 

NT1 = 2eNT -1 
IM2 = Mn2 

ANT] = NT1 

AIM2 = IM2 


SP(N) = SIDE*ANT1/((AIM2:4+2 )+Q) 
SP(M) = SP(M) + SIDE#(2.*AIM2-1. )/( (ATM2s+#2 )Q) 


CONTINUE 

DO 101 N = 3, 300 
MLESS1 = M-1 

DO 130 I =1,M 


DL(I) = DM(z) + VM(I)#12,*T 
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DO 131 I = 1,MLESS1 

CL(I) = DL(I) - DL(Z+1) 
IF(DL(M)) 12,12,1) 

be = 0.0 

GO TO 2) 

IF(DL(M) - Q) 16,16,18 

Ba = 0,0 

GO TO 2h 

VALUE = DL(M) - Q 

IF (VALUE-DE) 22,22,20 

DE = VALUE 

GO TO 2h 

DE = DE 

DP(N) = DE 

R(M) = (DL(M)-DP(N) )#*SP(M)s#(1.+ZVM(M) ) 
IF (DP(N)) 25,25, 26 

IF (DP(N) - DP(N-1)) 197,197, 25 
CONTINUE 

VAL = 8(1)#CL(1) 

IF (CL(1)-CM(1)) 5,5,h 
Pel) = VAL 

GO TO 33 

CL(1) = CM(1) 
FL(1)=VAL/(RES1i#2 )=(1, /(RESL#2 )-1. )*S(1)#CL(1) 
VAL2 = S(2)#CL(2) 
IF(CL(2)-CM(2)) 35,35, 3 
FL(2) = VAL2 

GO TO 37 

G2) = CM(2) 

FL(2) = VAL2/(RES2%#2 )~(1. /(RES2242 )-1. )#S(2)#CL(2) 
IF(CAP) 38, 36, 36 

LF(FL(2)) 39, 36, 36 

FL(2) = 0.0 

CONTINUE 

IF(FL(1)) 8,133,133 

FL(1) = 0.0 

Bowe2 i = 3.MLESS] 

FL(I) = CL(I)*S(I) 
Meare )) 213,211,210 

IF (FTEN1L-FL(2))211, 211, 21) 
FTENL = FL(2) 

EOero 211 

IF (FMAX1-FL(2))215,211, 211 
FMAX1 = FL(2) 

CONTINUE 

MC =M/2 +1 

IF (FL(MC ) )216, 217, 218 

IF (FTEN2-FL(MC ) )217, 217,219 
FTEN2 = FL(MC) 

GO TO 217 

TF (FMAX2-FL(MC ))220, 217,217 
FMAX2 = FL(MC) 
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CONTINUE 
IF(FL(IML) )221, 222,22 

IF (FTEN3-FL( IML) )222, 222,223 

FTEN3 = FL(IM) 

GO TO 222 

IF (FMAX3-FL(IML) )225 , 222,222 

FMAX3 = FL(IML) 

CONTINUE 

IF (RUB)L9,U9, bh 

DO 48 I = 3, MLESS1 

DPSL(L) = DPSM(T) 

CHANGE = ae = 0 

SUM = DL(I) + 

TF (DPSL(Z) -CHANGE) 5. 46,46 

DPSL(I) = CHANGE 

IF (DPSL(I )-SUM)48,48, 17 

DPSL(L) = SUM 

R(L) = pe - DPSL(I )}#SP(T +(1. +ZP#VMCT )) 
DO 50 I =1,M 

Vig) = vM(T) + (FL(I-1) - FL(I) = R(I))#£*32.2/w(T) 
CONTINUE 

DO 75 K =1,M 

or 1, 

IF (VL(K)) 75,75, 800 

CONTINUE 

GO TO 199 

CONTI NUE 

IF(VL(2) + CHECK)99,99,190 

IF (VL(M)-CHECK )102, 102. 19h 


DO 98 K = 1,M 
DS(K) = DM(K) 
DM(K) = DL(K) 
DL(K) = 0,0 
VS(K) = VM(K) 
VM(K) = VL(K) 
VL(K) = 0.0 
CS(K) = CM(K) 
Gen) = CL(K) 
CL(K) = 0.0 
FS(K) = FM(K) 
FM(K) = FL(K) 
FL(K) = 0.0 
DPSS(K) = DPSM(K) 
DPSM(K) = DPSL(K) 
DPSL(K) = 0.0 
CONTINUE 

GO TO 196 


WRITE OUTPUT TAPE 6,191,N 
FORMAT (6H VELOCITY OF PILE CAP EXCEEDED TWICE THE RAM VELOCITY W 


HEN N WAS 13) 
GO TO 196 
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19, WRITE OUTPUT TAPE 6,195 ,N 
19S O FORMAT (64H VELOCITY OF PILE TIP EXCEEDED TWICE THE RAM VELOCITY W 


ny7 
198 


1 
200 
196 
607 
606 


305 
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HEN N WAS 13) 

GO TO 196 

WRITE OUTPUT TAPE 6,198,N 

FORMAT (26H DP BECAME CONSTANT AT N = 13) 
GO TO 196 

WRITE OUTPUT TAPE 6,200,N 

FORMAT (52H ALL VL WERE SIMULTANEOUSLY NEGATIVE OR ZERO AT N = Ik) 
CONTINUE 

WRITE OUTPUT TAPE 6,607,RU 

FORMAT (27H ULTIMATE RESISTANCE (RU) = E15.h) 
WRITE OUTPUT TAPE 6,606, (SP(L),L=1,™) 
FORMAT (17H NEW VALUES OF SP/(8E15.h)) 
WRITE OUTPUT TAPE 6,305,DP(N) 

FORMAT (7H SET = E10.) 

BLOW = 1./DP(N) 

WRITE OUTPUT TAPE 6,805, BLOW 

FORMAT (66H 

BLOWS/IN= E10.3) 

FTENL = FTENL/AREAL 

FTEN2 = FTEN2/AREA2 

FTEN3 = FTEN3/AREA3 

WRITE OUTPUT TAPE 6, 307,FTENL,FTEN2 ,FTEN3 


O FORMAT(30H MAXIMUM TENSION AT HEAD = E15.4,/ 
1 30H MAXIMUM TENSION AT CENTER = E15,)/ 
2 30H MAXIMUM TENSION AT TIP = E154) 


FMAX1 = FMAX1/AREAL 

FMAX2 = FMAX2/AREA2 

FMAX3 = FMAX3/AREA3 

WRITE OUTPUT TAPE 6, 306, FMAX1, FMAX2, FMAX3 
FORMAT (26H MAXIMUM FORCE AT HEAD = E15, L/ 
26H MAXIMUM FORCE AT CENTER = E15.h/ 

26H MAXIMUM FORCE AT TIP = £15.)//) 
CONTINUE 

RU = RU + 0000. 

CONTINUE 

IF (ITURN )205 , 205 , 207 

CAlMenpenT 

END 

DATA 
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APPENDTX. D 


DETAILS OF THE "HELEY FORMULA COMPUTER PROGRAM" 


‘General: This program is written in Fortran Language for the 
IBM 7090 computer, but may be used equally well with the 709 or 709 
computer. It has as its objective the systematic production of data of 
set vs. ultimate ground resistance using the Hiley type formula. The 
program will take given pile and hammer information and beginning with 
an ultimate ground resistance of 5 tons will range in increments of 
5 tons until an ultimate ground resistance of 2:0 tons is reached. At 
each value of ground resistance the set is computed and written out. 
The program 1s arranged so that any type of hammer or pile can be used. 
the coefficients for use with the Hiley formula are as recommended by 
Chellis (ref. 5) and are built into the program except as indicated 
under input data preparation section below. 

Termination: No conditions for termination are built intc 
this program. If extensive use were planned for this program, an "if" 
statement should be added so that it will cause the calculations to 
terminate when the computed set becomes negative. 

Input Data Preparation: All data is punched on standard data 
cards in the exact order indicated below. Upon completion of punching 


the data deck, it is placed immediately after the « DATA card of the 


source program, 


CfRD ORDER DESCRIPTION 
rE If another set of data follows this set for 


another case, punch +1 in spaces 9 and 10; if 
not, punch -1. Use no decimals. 


‘ — «s 
el a lee “a -_— > 9S 







Eee ) tte (2 oye as 
= Ce | ae eye me pa Cee Ce 
=" - eye a Oe cee -ae aad ti 








ee! ee - =—_ ate Sore omy I we 
— — er eee See ee ali penne eet Ce = 









» | — lie . weiss >, 
- s A 4M re « hie » Tee 
i ee ee em ——_— oo oe 


-—| a -_ oo eee) «) Gee 


“s —— Pt we oer Ta 
Se >» —_— © oe we wee fe . pel 
—_ * 2[2————— |) a fs en 
; | aI) ge Hh &§ eae 
- » « meee (> ee « 
. ‘= \ ~— am @ C84 one ee 
‘i <4 


nee _ a 


- te Ne AS PT ee ee 
-_—i i= ~T & t a. -—s 
- 78 ‘ = ' . im, - 
=a ¥ 
GED, OOS 
—_ ‘on 4s ° ; 


LO 


deli 


ie 


aS 


1h 


146 


if the hammer is double acting, punch +1 in 
Spaces 7 and 10; if not; punch -1. Use no 
decimals. 


if the pile is steel, punch +1 in spaces 9 and 
10; if wood, -1 in spaces 9 and 10; if concrete, 
punch O in space 10. Use no decimals, 


If triangular side distribution is desired, punch 
+] in spaces 9 and 10; if rectangular side dis- 
tribution is desired, purch -1 in spaces 9 and 
nO, 


in the first 10 spaces punch the hammer effi- 
ciently using a decimal point. 


In the first 10 spaces punch the weight of the 
ram in pounds using a decimal point. 


In the first 10 spaces punch the coefficient 
of restitution using a decimal point. 


In the first 10 spaces punch the weight of the 
pile including the shoe and driving cap for 

drop hanmers and single acting steam hammers. 
Weight of pile including weight of anvil in case 
of double acting or differential acting steam 
hammers, All in pounds and with a decimal point. 


In the first 10 spaces punch the cross sectional 
area in square inches of the pile at its head 

or the area of the capblock if one is used. Use 
a decimal point. 


same as card 9 but for the cross sectional pile 
area at mid-length of the pile. 


Same as card 9 but for the cross sectional pile 
area at the pile tip. 


Punch in the first 10 spaces the ground quake 
in inches, using a decimal point. 


Punch in the first 10 spaces the modulus of 
elasticity in pounds per square inch using a 
decimal point. 


Punch in the first 10 spaces the pile length in 
feet using a decimal point. 


lll 
‘ 
’ ee 
=F 
' 
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a Punch in the first 10 spaces the height of free 
fall of the ram, in inches, for drop hammers 
or 


the normal (shortest) stroke of ram in inches 
for single acting steam hammers 
for double acting, Pre eee acting steam and 
diesel hammers use the rated energy per blow in 
foot pounds as published by the manufacturers. 
Running the Program: The data deck of cards just described is 
placed immediately after the source program deck and is ready for 
delivery to a machine operator for processing. 
Results: In the output will be repeated the original input 
information, values of the ultimate ground resistance, per cent of 
joad carried by the point, and the set. The results are printed in 
decimal form with a multiplier so that 0,6192 EOO, for example, means 
0.6192 x 10°, or 0.6192. 
The Program: The complete program is presented so thab any who 
wish to use it may do so. Each line of the program should be punched 
on a separate card and assembled in the exact order indicated. Anyone 


preparing cards from this print-out is cautioned to verify the program 


with test cases before relying on the results obtained. 
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D REESE C2165 CIVIL ENGR’ G HILEY 7 150 
XEQ 
NOBIN 
LABEL 
HILEY FORMULA PROGRAM 


206 READ INPUT TAPE 5,20, ITURN, IEQUIP,IMAT, MUD 


20 FORMAT(I10) 
READ INPUT TAPE 5,21,E,WR,RES,WP,APH,APAV, APT,Q,ELAS,PL,H 
21 FORMAT (F10.)) 
IF(IMAT) 1,3,5 
1 WRITE OUTPUT TAPE 6,2 
2 FORMAT (12H WOODEN PILE) 
GO T0 7 
3 WRITE OUTPUT TAPE 6,), 
)} FORMAT (14H CONCRETE PILE) 
GO TO 7 
5 WRITE OUTPUT TAPE 6,6 
6 FORMAT(11H STEEL PILE) 
7 IF(IEQUIP) 8,12,10 
8 WRITE OUTPUT TAPE 6,9 
9 FORMAT(3LH SINGLE ACTING OR DROP HAMMER USED ) 
GO TO 12 
10 WRITE OUTPUT TAPE 6,11 
11 FORMAT(26H DOUBLE ACTING HAMMER USED) 
12 IF(MUD) 13,205,15 
13 WRITE OUTPUT TAPE 6,14 
1), FORMAT(33H SIDE DISTRIBUTION IS RECTANGULAR ) 
GO TO 17 
15 WRITE OUTPUT TAPE 6,16 
16 FORMAT (32H SIDE DISTRIBUTION IS TRIANGULAR ) 
17 WRITE OUTPUT TAPE 6,18,E 
18 FORMAT(13H EFFICIENCY= E10, h) 
WRITE OUTPUT TAPE 6,19,WR 
19 FORMAT (16H WEIGHT OF RAM= E10, h) 
WRITE OUTPUT TAPE 6,22,RES 
22 FORMAT (29H COEFFICIENT OF RESTITUTION= E10, hh) 
WRITE OUTPUT TAPE 6,23,WP 
23 FORMAT (13H WT OF PILE= E10, Ih) 
WRITE OUTPUT TAPE 6,2h, APH 
2h FORMAT (26H AREA PILE HEAD IN SQ IN= E10. 4) 
WRITE OUTPUT TAPE 6,25, APAV 
25 FORMAT (20H AVERAGE PILE AREA= E10.) 
WRITE OUTPUT TAPE 6,26, APT 
26 FORMAT (16H AREA PILE TIP= E10,)) 
WRITE OUTPUT TAPE 6,27,Q 
27 FORMAT (8H QUAKE= E10, fh) 
WRITE OUTPUT TAPE 6,28,ELAS 
28 FORMAT (17H YOUNGS MODULUS= F10,);) 
WRITE OUTPUT TAPE 6,29,PL 
29 FORMAT (20H PILE LENGTH IN FT= E10.)) 
IF (IEQUIP) 31,35, 33 
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31 WRITE OUTPUT TAPE 6,32,H 
32 FORMAT (24H HAMMER DROP IN INCHES= E10, h) 
GO TO 35 

33 WRITE OUTPUT TAPE 6, 3h, H 
3) FORMAT (OH RATED ENERGY OF HAMMER IN FOOT POUNDS= E10. hh) 
35 RU = 10000, 

DO 20h J = 1,8 

WRITE OUTPUT TAPE 6, 30,RU 
30 FORMAT (20H THE VALUE OF RU IS E10.h//) 


Pl = RU/APH 
P2 = RU/APAV 
P3 = _. 


DO 203 I = 1,3 
IF(I-1) 60, 0, 2 
hO PLMOD = PL 
WPC = WP/2, 
WRITE OUTPUT TAPE 6,1 
1 FORMAT (29H ALL OF LOAD CARRIED BY POINT) 
GO TO 60 
42 IF(I-2) 60,4452 
hl WRITE OUTPUT TAPE 6, 7 
h? FORMAT(31H ONE HALF LOAD CARRIED BY POINT) 
IF(MUD) 46,60,18 
6 PLMOD = O, 75 xPL 
WeC = WP 
GO TO 60 
48 PLMOD = h, « PL / 6, 
WPC = WP 
GO TO 60 
52 IF(I-3) 60, 54, 60 
Si, WRITE OUTPUT TAPE Onaoe 
SS FORMAT(29H ALL LOAD CARRIED BY FRICTION) 
IF(MUD) 56,60,58 
56 PLMOD = PL/2 
WPC = WP 
GO TO 60 
58 PLMOD = 2, * PL / 3. 
WPC = WP 
60 C3 = Q 
IF(IMAT) 70, 80,90 
70 Cl= 0.00010 * Pl 
C2 = 0,00008%P2%*PLMOD*15 00000. /ELAS 
GO TO 100 
80 Cl= 0,00025 * PL 
C2 = 0,00000))+P2%PLMOD#+3000000. /ELAS 
GO TO 100 
90 Cl = 0.00008 * Pl 
C2 = 0,000000h+P2*PLMOD*30000000, /ELAS 
100 DP=(( (ExWR3H )+(WR+WPCxRES#:2 ) )/( (WR+WPC )*RU) )-(C1+C2+C3)/2. 
IF (IEQUIP) 104,203,102 
102 DP = DP + (c14024+03)/2. 
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DP 12.%DP 

DP = DP ~(C1+02+C3)/2. 
104 BLOW = 1./DP 

WRITE OUTPUT TAPE 6,201,DP 
201 FORMAT(8H SET IS E20.h) 

WRITE OUTPUT TAPE 6,202, BLOW 
202 FORMAT(20H BLOWS PER INCH ARE E20,)/) 
203 CONTINUE 
20) RU = RU + 10000, 

IF (ITURN) 205,205,206 
205 CALL EXIT 

END 

DATA 





Vea 


APPENDIX & 


This appendix includes some of the graphs of computer solutions 
which were used both for correlation with actual pile cases and for 
studying the effect of varying the size of hammer, length of pile, 
ground quake, point damping, and side damping. On each sheet several 
curves are drawn for various values of frictional resistance expressed 
as a percentage of the total ultimate bearing capacity of the pile and 
indicated as follows: 

Curve A -- End bearing; no side friction 

Curve B ~~ 75 per cent end bearing; 25 per cent side friction. 

Curve C -- 50 per cent end bearing; 50 per cent side friction. 

Curve D -- 25 per cent end bearing; 75 per cent side friction. 

Curve E -- 100% friction. 

Rectangular side distribution was used. The description of the driving 
equipment, pile, and soil properties for which the computations were 
made appear on each sheet, but for convenience of use in future correla- 
tion work an index is included which summarizes this information for 
each sheet. ‘The figure numbers and page numbers are identical in this 


appendix. 
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| INDEX OF FIGURES IN APPENDIX E 
1. Piles Driven by a Vulcan #1 Hammer 
a. O-foot 12BP53 Piles: 
Figure 

Point Damping Side Damping Ground Quake Number 
O20 0.0 OFL0 Rl 
Q.10 0,033 0.10 K2 
0.15 ORC O.10 E3 
0.20 0.067 OL Ely 

0.30 0.10 0.10 ES 

| 0.40 O33 Orie E6 
G50 0.167 Ooo E7 
CLO 0.033 0.20 E8 
0.20 0,067 0,20 E9 
0.30 QO 10 Oreo E1O 
0.40 Ons. 0.20 E11 
0.50 0,167 0220 E12 
@e20 0,033 0.30 E13 
OF20 0.067 Oe El 


b. 80-foot 12BP53 Piles: 


r 0.0 0.0 0.10 E15 
0.10 0,033 OMe E16 
0.20 0.067 0,10 E17 
0.30 0.10 0.10 E18 
0.40 0.133 O10 E19 


0.50 0.167 0.10 B20 


=~... - i 
— tiie 





(80-foot 12BP53 Piles cont'd. ) 


C. 


Point Damping 


0,10 
0.20 
0.30 
0.40 
0.50 
0.10 


0.20 


0.0 

0,10 
0.20 
0,30 
0.40 
6250 
0.10 
0.20 
O20 
0.40 
0.50 
0.10 


OF 20 


0.15 


Side Damping 


0,033 
0.067 
0.10 

O51) 
Onloy, 
0,033 
0.067 


120-foot 12BP53 Piles 


0.0 
0,033 
0.067 
0.10 
Orr33 
0.167 
0.033 
0.067 
0.10 
Orkes 
Onley 
0.033 
0.067 


0.0 


Ground Quake 


O20 
0.20 
O20 
0.20 
Oncd 
0550 


0.30 


OL ue 
O720 
0.10 
Oko 
0.10 
0,10 
Orzo 
0,20 
0.20 
0.20 
0,20 
0,30 


0.30 


18-inch~square concrete pile 0 feet long: 


0,10 


Figure 


Number 


EL 
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2. Piles Driven with a Vulcan O Hammer 


12BP53 Piles 80 feet long: 


Figure 
Point Damping Side Damping Ground Quake Number 
720 0,067 0.10 Eh2 
0.4.0 Ooes 0.10 Eh,3 
12BP53 Piles5 feet long: 
0.10 07033 Oso Eu 
1)4BP89 Piles 5 feet long: 
0,15 0.05 0,10 ES 
6.20 0,067 0.10 E6 
0.20 0.067 0.20 Eh? 
1iBP117 Piles lS feet long: 
O15 B05 O70 ELS 


Oac0 0.067 Oec0 EY 
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